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The physics of graphene is acting as a bridge between quantum field theory and condensed matter 
physics due to the special quality of the graphene quasiparticles behaving as massless two dimen- 
sional Dirac fermions. Moreover, the particular structure of the 2D crystal lattice sets the arena to 
study and unify concepts from elasticity, topology and cosmology. In this paper we analyze these 
connections combining a pedagogical, intuitive approach with a more rigorous formalism when re- 
quired. 
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I. INTRODUCTION 

One of the sad consequences of the very fast development of physics during the last half of a century is that it 
has started to loose its unity. Especially, the gap between condensed matter physics which is a key ingredient of our 
scientific understanding a "world around us" and microphysics investigating the most fundamental laws of nature 
(high-energy physics, quantum field theory, gravity and cosmology, etc.) has become wider. That is why examples of 
fruitful exchange of ideas and methods between these two branches of physics are very important now. A prototype 
example of the relevance of macrophysics for microphysics is the idea of spontaneously broken symmetry which first 
appeared in Landau's theory of second-order phase transitions [1 and then turned out to be one of the most important 
concepts in both condensed matter [2 and fundamental physics [3, 4 . Conversely, the idea of the renormalization 
group, being developed first to study the problem of infinities in quantum electrodynamics (QED) has been of crucial 
importance to solve problems of condensed matter physics as difficult as critical behavior and the Kondo problem [5^9^ . 
The use of Mossbauer effect to check the general relativity theory [10] gives an amazing example of an entanglement 
of condensed matter physics, nuclear physics, and gravity. 

One of the fields where condensed matter physics meets quantum field theory and cosmology is that of the super- 
fiuidity of helium 3 [11 . A recent development in material science provides a new and unexpected bridge between 
condensed matter and high-energy physics. The experimental discovery of graphene, a two-dimensional allotrope of 
carbon formed by a single carbon atom sheet, was made in 2004 when a technique called micromechanical cleavage 
was employed to obtain the first graphene crystals [HI [13] . The observation of a peculiar "Dirac-like" spectrum of 
charge carriers and an anomalous quantum Hall effect in graphene [TH [15] has ignited an enormously growing interest 
to this field (for a review, see [16 - 20 ). One of the most interesting aspects of the graphene physics from a theoretical 
point of view is the deep and fruitful relation that it has with QED and other quantum field theory ideas |2Tti29] . 
Probably the most clear example is the Klein paradox [30l [31] , a property of relativistic quantum particles of being 
able to penetrate with a probability of the order of unity through very high and broad potential barriers. Previously it 
was discussed only for experimentally unattainable (or very hard to reach) situations such as particle-antiparticle pair 
creation at the black hole evaporation, or vacuum breakdown at collisions of super-heavy nuclei (a rather complete 
reference list can be found in [32 ). At the same time, it appeared to be relevant for graphene-based electronics and 
for electronic transport in graphene [25l [28] [33H35I . 

Here we focus on a particular aspect of the physics of graphene that establishes its connection with QED, field 
theory and gravity, namely, the appearance of gauge fields in graphene and their effect on its properties. The concept 
of a gauge field in general has been extensively discussed in condensed matter physics, especially, in relation with 
modeling different types of topological defects, phase transitions, and properties of glasses [36-^40] . In the modern 
context of graphene the gauge fields were first introduced in (HJ [42] in relation with the problem of weak localization. 
It is well-established by now, both experimentally and theoretically [4T1 [I3ti47] , that graphene is always corrugated 
and covered by ripples which can be either intrinsic [38| |43l [441 HZ] or induced by a roughness of substrate [45l [46] . In 
general, the departure from fiatness of graphene leads to the appearance of an inhomogeneous pseudomagnetic gauge 
field [m [42] acting on the charge carriers. This has important consequences affecting the character of quantum Hall 
effect in graphene, it produces an effective source of charge carrier scattering, and it can provide a mechanism for an 
observed charge inhomogeneity, as discussed in detail below. Whereas a smooth deformation of the graphene sheets 
produces a gauge field similar to the electromagnetic one, different topological defects in graphene inducing inter- 
valley (Umklapp) processes can be considered sources of a non- Abelian gauge field [48] [49] . On the other hand the 
description of some topological defects as cosmic strings [50] creates an interesting connection with general relativity. 
This will be also described in the present review. 

The work is organized as follows. Section |ll| presents the electronic features of the model that are most relevant 



for the physics discussed. Section |IIA| gives a short description of the way in which the Dirac equation approximates 



the electronic states of graphene in the long wavelength limit emphasizing the need for two component wavefunctions. 



In subsection II B we show the topological stability of the Fermi points in the honeycomb lattice towards small 



perturbations like lattice deformations or disorder. Then, in subsections II C and |IID| we describe the way in which 
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FIG. 1: (Color online) Lattice and unit cell for monolayer graphene in real space. 



gauge fields were defined in physics and in geometry and we show the different types of effective gauge fields that 
appear in graphene associated to various disorder types. Different physical mechanisms which can contribute to these 
gauge fields are discussed next, in sections III (topological defects) and |IV| (smooth deformations). Observable effects 
related to the gauge fields are described in sections [V| (microscale gauge fields) and [Vl] (macroscale gauge fields) . 



Section VII presents a summary of the work, open questions, and some possible future developments. Technical 



aspects related to general features of spinors in a curved two dimensional surface are explained in Appendix [A] and 
some auxiliary material from the elasticity theory is presented in Appendix [B] 



II. THE HONEYCOMB LATTICE: SPINORS AND GEOMETRY IN TWO DIMENSIONS 

A. The low energy electronic excitations of graphene. Continuum model 

Monolayer graphite - graphene - consists of a planar honeycomb lattice of carbon atoms shown in Fig. [l] In the 
graphene structure the in-plane a bonds are formed from 2s, 2px and 2py orbit als hybridized in a sp^ configuration, 
while the 2pz orbital, perpendicular to the layer, builds up covalent bonds, similar to the ones in the benzene molecule. 
The (J bonds give rigidity to the structure, while the tt bonds give rise to the valence and conduction bands. The 
electronic properties around the Fermi energy of a graphene sheet can be described by a tight binding model with 
only one orbital per atom, the so-called 7r-electron approximation because there is no significant mixing between 
states belonging to a and tt bands in 2D graphite. Within this approximation a basis set is provided by the Bloch 
functions made up of the 2pz orbitals from the two inequivalent carbon atoms A and B which form the unit cell of the 
honeycomb lattice. One may define two Bloch wave functions to be used in a variational (tight-binding) computation 
of the spectrum 

^i{K) = ^e'^-^''^''^^{f-fi-^ , i = A,B (1) 
t 

where the sum runs over all the points in the direct lattice, i.e., t = 711^1+712^2, (r^, vb) are the positions of the atoms in 
the unit cell, and ^{f) is a real (yr-type) atomic orbital. As it is well known, a simple tight-binding computation [5T| [52] 
yields at the neutrality point a Fermi surface made of the six single points located at the corners of the hexagonal 
Brillouin zone. Only two of them are inequivalent and can be chosen to be located at Ki = —2gi/3 — g2/3 and 
K2 = —Ki, where ti • gj = 27r6ij. A low energy expansion around any of the two Fermi points gives an effective 

hamiltonian linear in momentum which reduces to the massless Dirac equation in two dimensions derived from the 
Hamiltonian: 

no = -hvF J d^r^{r){ia:,da: + iaydy)^{r) , (2) 

where cr^ are the Pauli matrices, = (3ta)/2, t is the nearest-neighbor hopping parameter, and a = 1.42 A is the 
distance between nearest carbon atoms. 

The components of the two-dimensional spinor: 

^ ( ) (3) 
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correspond to the amplitude of the wavefunction in each of the two sublattices which build up the honeycomb structure. 
The two Fermi points of the honeycomb lattice give rise to two similar equations related by time reversal symmetry. 

We can combine the two spinors attached to each Fermi point into a four component Dirac spinor and write the 
four-dimensional Hamiltonian: 

Hd = -ivFh{l (g) (Jidx + Ts (g) cr2dy), (4) 

where a and r matrices are Pauli matrices acting on the sublattice and valley degree of freedom respectively. 

Summarizing, the internal degrees of freedom of charge carriers in graphene are: Sublattice index (pseudospin), 
valley index (flavor) and real spin, each taking two values. The real spin is irrelevant for the issues discussed in this 
review and will not be taken into account except for an additional degeneracy factor 2 in some quantities. 



B. Topological stability of the Fermi points 

The properties of graphene presented in this work as well as most of the exotic properties of the material lie on 
the special character of its low energy excitations obeying a massless two dimensional Dirac equation. Hence it is 
important to establish the robustness of the low energy description under small lattice deformations, disorder and 
other possible low energy perturbations. In this section we will analyze the existence and stability of the Fermi points 
in the graphene lattice. A related issue is the possibility to open and control a gap in the material, crucial for electronic 
applications. We will analyze under which conditions a gap can open in neutral graphene. We will see that the Dirac 
points in graphene are topologically preserved under rather general circumstances and Coulomb interactions, stress, 
phonons, and other moderate perturbations of the lattice will not open a gap. The analysis of this section follows 
closely ref. [53 . 

The discrete symmetries of the system playing a very important role in the analysis we will summarize them first. 
As it was described in Section II A[ the Fermi surface of neutral graphene consists of two independent Fermi points 



that can be chosen in the Brillouin zone to be located at K and -K. The effective low energy Hamiltonians around 
these points are 

/ Q ^* \ 

H{Ki + ^) ^ ( ^ Q ) = kxCTx + kyay, (5) 

and 

^ ^ / —k \ 
H{-Ki ^k) ^ i Q j -/ca^o-^ + kyCFy. (6) 

where k = kx -\- iky and are the Pauli matrices. 

A gap will be opened by a /c-independent, translationally invariant perturbation of the type 

where a = ax -\- iay. The spectrum becomes E = zb^a^ -\- \k -\- ap and a gap 2|a2|is generated. However as we will 
see, such a perturbation is not allowed if the discrete symmetries of the system are to be respected. The two relevant 
discrete symmetries in graphene are time-reversal T : t ^ —t and spatial inversion / : {x^y) {—x^ ~y)- The reality 
of the TT orbitals implies that time reversal merely reverses K 

T^,{K) = = (8) 

whereas the spatial inversion also exchanges the two types of atoms 

I^a{K) = ^b{-K) , I<^b{K) = ^a{-K)- (9) 
Invariance under these symmetries imposes the following constraints on the Hamiltonians: 

T: H{K) = H*{-K) 

I : H{K) = (JxH{-K)(Jx. (10) 
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From eqs. (I9| and (10) we see that none of the two symmetries are obeyed for the effective description around a single 



Fermi point. The graphene system is only invariant under T and / individually if both Fermi points are considered 
simultaneously. Nevertheless the product TI imposes a constraint on the form of H{K) at each given Fermi point: 



TI: H{K) = a^H*{K)a,. 



(11) 



If the allowed perturbation respects TI this implies that IIii{K) = i^22(^), what enforces =0 in ([7|) and no gap 
opens. 

This has an interesting topological interpretation, which extends the previous arguments to /c-dependent — but 
translationally invariant — perturbations. The low energy hamiltonian II{Ki + ^) in (jsl) defines a map from the circle 



ky = R to the space of 2 x 2 Hamiltonians H = h • a: 



k = i?e'^ {hx,hy,hz) = i?(cos6>,sin6>,0). 



(12) 



Since Fermi points correspond to zeroes of the determinant —Det{II) = h^. -\- hy -\- h^, a perturbation will be able 
to create a gap only if the loop represented by the map (12) is contract ible in the space of Hamiltonians with non- 
vanishing determinants, which is just — {0}. This is clearly the case, since 7Ti{R^ — {0}) 
other hand, Hamiltonians invariant under TI are represented by points in i?^, and we have 



7ri(6'2) = 0. On the 



7r,{R'-{0})=7r,{S') = Z 



(13) 



This means that non-trivial maps such as the ones implied by ([5| can only be extended to the interior of the circle 
by going through the origin, i.e., by having at least one zero. This precludes the creation of a gap. The non-trivial 
topological charge associated to the topological stability of each individual Fermi point in graphene is simply the 
winding number which can be computed from 



47rz Jq 



(14) 



It is important to realize that the two Fermi points have opposite winding numbers implying that they could annihilate 
mutually if brought together by a perturbation. 

The nontrivial topological charge associated to the topological stability of individual Fermi point in very general 
2D systems has been discussed in [54 and in [ 55H57] in the context of topological insulators. 



C. Gauge fields in physics and geometry 

Gauge invariance plays a key role in the quantum field theory (QFT) description of fundamental forces between 
elementary particles. The seemingly abstract concept of a non-Abelian gauge field introduced first in QFT to describe 
the electroweak interaction followed by the experimental discovery of the W and Z bosons is one of the most impressive 
achievements of theoretical physics. Before introducing the various gauge fields appearing associated to the physics 
of graphene and in order to clarify their specific nature we will make a brief description of the classical concept of 
gauge invariance and of the associated gauge fields. 

The concept of gauge invariance emerged from classical electrodynamics. In particular the electromagnetic field 
(E,B) is expressed in terms of potentials (<l>. A) through: 

E = -(V<I> + S^A) , B = VAA. (15) 

The fields do not change under the transformation 

A^A + Vx , ^^^-dtx. (16) 

where x is an arbitrary (smooth) function of space. This invariance was shown to survive to the Quantum Mechanics 
description of a charged spinless particle in an electromagnetic field provided that the wave function was simultaneously 
transformed to 

^^^exp(iex). (17) 

The relativistic wave equation for a spinless particle with charge e interacting with electromagnetic fields is derived by 
first performing the substitution ^ — eA^, where = {A^ = A) is the 4- vector electromagnetic potential 
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and then performing the usual substitution ihd^. A formal solution for the wave function of a particle interacting 

with the electromagnetic potential can be written in terms of the solution without interaction as 



^ = exp 



ie / A^dx. 



^0. 



(18) 



The original discovery by V. Fock was that the quantum dynamics, i. e. the form of the quantum equation, remains 
unchanged by the transformations (16) if the wave function of the particle is multiplied by a local (space-time- 
dependent) phase. The modern description and the term " gauge invariance" was established by H. Weyl. An 
interesting account of the historical development of the gauge concept is given in ^ . 

The first example of a QFT gauge model is four dimensional QED. A free spin 1/2 Dirac fermion with charge e 
and mass m is described by the action 



d'^x'i [7^5^ + m] 



which is invariant under the global f/(l) group transformations: 



*(a;) U'^{x) , *(a;) -> f/*\E'(a;) , U = exp(iex), 



(19) 



(20) 



where x is a constant. Gauge invariance requires invariance of the action under the local group of transformations 
obtained by replacing x x(^)- This can be achieved by replacing the derivative in (19) by the covariant derivative 



Djj, = djj, -\- ieAjj^: Under a local U{1) transformation defined by (20) with a space-time dependent function x{^)^ 
A^{x) transforms as A^ A^ — d^x-> ^ generalization of ([16|. 

The invariance of the Maxwell equations under special relativity allows a formulation of (classical) electromag- 
netism in terms of quadrivectors and tensors. The equations can be written in a covariant way by introducing the 
electromagnetic tensor F^y defined by 



^ij k -^k 1 



(21) 



and the quadricurrent = (p, J) made of the charge density and the current. In terms of these geometric objects 
the four Maxwell equations reduce to 



dxFa 



0, 



(22) 



S^F^^ = (23) 

The conservation of the current dj^J^ =0 follows from the antisymmetry of F^j^. The first equation is identified as a 
Bianchi identity that can be integrated by introducing a gauge field A^ such that F^jy = d^Aj^ — djyA^. It is readily 
verified that two gauge fields related by the gauge transformation A^^ = Aj^ — dj^Q give rise to the same electromagnetic 
tensor field. Maxwell's equations can be derived from the action 

S{A,J) = jd^xlF^^F'^'^ + Jf^A''], (24) 

which is precisely the full QED action. 

The concepts of gauge fields and covariant derivatives can be translated into the language of differential geometry, 
based on differential forms. In general, the gauge field has a mathematical interpretation as a Lie-valued connection 
and is used to construct covariant derivatives acting on fields, whose form depends on the representation of the group 
under which the field transforms (for global transformations) . The field tensor F^^^ is a curvature 2-form given by the 
commutator of two covariant derivatives. It is an element of the Lie algebra associated to the gauge group. The gauge 
connection generates parallel transport of the geometric objects under gauge transformations. The generalization of 
U{1) to non-Abelian groups as SU{N) is straightforward, the main modification arises in the definition of the field 



strength (21) that becomes F^j, = — d^A^ + [A^, A^^]. Since most of the gauge fields arising in the graphene 

context will not have dynamics we will not discuss this point further. 

Einstein's General relativity can be also interpreted as a gauge theory where gauge invariance is invariance under 
diffeomorphisms (local smooth changes of coordinates) in the space-time manifold. The connection which generates 
parallel transport plays the role of the gauge field. The similarity is better appreciated in the vielbein formalism 



introduced in Section III C and Appendix 1 where a gauge transformation corresponds to a change of local frame (a 
local Lorentz transformation). Gauge invariance corresponds to the independence of field equations from the choice 
of the local frame. The spin connection plays the role of the gauge field. 



7 



In any gauge theory physical observables are related to gauge- invariant operators. The gauge invariance allows to fix 
some conditions on the gauge potentials that will not affect the physical properties. In quantum gauge theories gauge- 
fixing is a delicate issue that will no concern us here. In classical electromagnetism, the gauge-fixing problem is simply 
the problem of choosing a representative in the class of equivalent potentials, convenient for practical calculations or 
most suited to physical intuition. In non-relativistic problems one of the most popular choices is the Coulomb gauge: 
V.A(t,x) = whose relativistic counterpart 9^A^(t,x) = 0,/i = (0,1,2,3) is the Landau or Lorentz's gauge. The 
freedom to choose a gauge condition is related with the full gauge invariance of the action. When fictitious gauge 
fields are generated by analogy with the gauge formalism but there is no dynamics associated to them it can happen 
that the gauge potentials are fixed by the physics involved and no extra conditions can be imposed. A particular 
example is provided by the strain fields discussed in Sect. |IV A[ Gauge fields were introduced in condensed matter in 
the early works of refs. [59l [60]. 



D. Types of gauge fields 



Different physical mechanisms give rise to perturbations in the Dirac equation described above which can be 
described mathematically as effective gauge fields. By definition, a gauge field appears in the off diagonal elements of 
the Dirac equation, which correspond to the hopping between sublattices. There are two ways in which this change 
in the hopping between sublattices can take place: 

• The topology of the lattice itself can change, requiring a redefinition of the two sublattices. This happens 
in the presence of topological defects, such as pentagons and heptagons (disclinations), or pentagon- heptagon 
pairs (dislocations). If the carbon atoms remain with threefold coordination, the local electronic structure is 
not altered very much, and an effective Dirac equation can be defined locally. The defect changes the global 
properties of the lattice. Modifications of the internal degrees of freedom of the electrons when they propagate 
over long distances can be described by gauge fields, as defined more precisely in the following section. 

The substitution of hexagons by rings with fewer or more sites, when the bond lengths remain more or less 
unchanged, leads to the curvature of the structure, according to Euler's theorem. As described below, these 
defects induce, in the continuum limit, the spin connection which modifies the Dirac equation of a spinorial field, 
extensively discussed in quantum field theory in curved spaces. Because these gauge fields arise from general 
topological features of the system, these gauge fields do not depend on any material parameter. 

• The hopping between tt orbitals in different sublattices can also change because the distance between the orbitals 
is modified, or because a change of symmetry allows for indirect hoppings through the a orbitals [61 , also present 
in the system. These fields depend on microscopic details of the material. The same changes in the electronic 
spectrum are induced when the atoms oscillate around their equilibrium positions, so that the parameters which 
describe the strength of these gauge fields contribute also to the electron-phonon coupling. 

Small static deformations of the lattice are described by the strain tensor. A field with the symmetries of a 
vector can be obtained by contracting the strain tensor with a third rank tensor. The trigonal symmetry of the 
honeycomb lattice allows for such a third rank tensor [62 , 63 (note that a vector gauge field cannot be defined in 
a fully isotropic system). The curvature away from the flat configuration can be described in a similar way by 
the curvature tensor. The induced gauge field in the Dirac equation is obtained in analogous way by contracting 
the curvature tensor and the third rank trigonal tensor described above. 



III. TOPOLOGICAL DEFECTS IN GRAPHENE 



As described in Sect. II A graphene consists of a planar honeycomb lattice of carbon atoms determined by the sp^ 
hybridization of the 2s, 2px and 2py orbitals of the carbon system. These a bonds give rigidity to the structure, while 
the TT bonds give rise to the valence and conduction bands. The elastic properties of the material are related with the 
a bonds and involve energies of the order of 7-10 eV similar to the bandwidth of graphene (~ lAeV). The carbon 
bond in graphite is one of the strongest chemical bonds occurring in nature. The Young's modulus of graphene is of 
the order of 350 N- m~^, one of the highest values known for any material j64]466] . 

Topological defects in solids are usually described by complicated boundary conditions in elasticity theory. It was 
very early realized that they may be described more simply as sources of a gravity-like deformation [40 . In this 
approach the boundary conditions imposed by defects in elastic media are accounted for by introducing curvature 
in the given space. In the continuum limit the crystal is described by a manifold where curvature and torsion 
are associated to disclinations and dislocations in the medium respectively. The Burgers vector of a dislocation is 



FIG. 2: (Color online) Honeycomb lattice of monolayer graphene. 



associated to torsion, and the Frank angle of a disclination to curvature. The combined effect of curvature and 
elasticity has been studied in [67 and the scattering of phonons in the presence of a disclination and of multiple 
disclinations in [68]. The fact that elastic deformations of a crystal lattice induce gauge fields that couple to the 
electronic degrees of freedom has been known for a long time [36] and it has been applied to nanotubes and graphene 
in [4r, 42 , 68^. The mechanics of defects of various classes in graphene has been analyzed recently in \69 . The authors 
discuss the failure of classical elasticity to describe divacancies and Stone- Wales defects. The defect formation under 
strain was studied in the context of carbon nanotubes in [70] . 

A. Topological defects. Formation and naturalness 

A very natural way of producing local curvature in the graphene lattice is by substituting some of the hexagons 
by pentagons (positive curvature) or heptagons (negative curvature). These types of topological defects are present 
in all the previously existing graphene structures (fullerenes and nanotubes), have been observed experimentally [711 - 
[73] and their elastic and electronic properties have been studied at length [74 . It is well known for instance from 
studies of stability of carbon nanotubes that vacancies and similar defects produced in the graphene structure by ion 
irradiation of the samples are mended by forming higher membered rings. Subsequently the unstable high-membered 
rings have been found to disappear by Stone- Wales transformations [75] , thus leading to very stable structures mainly 
constituted of five-, six-, and seven- membered rings [76^. High-resolution Transmission Electron Microscopy (TEM) 
with atomic sensitivity reported the first direct imaging of pentagon-heptagon pair defects in a single wall carbon 
nanotube heated at 2,273 K in [77] . 

The recent TEM experiments performed in suspended graphene sheets [43j [44^ seem to indicate that the observed 
structure can not be explained by strain-free deformations of graphene. Topological defects can be produced in the 
process of mechanical cleavage which involves quite high energies. They would induce strain- free curvature compatible 
with the experimental data. Although the observed tilt reported is argued to produce a strain too small to give rise 
to nucleations of defects [43], it can well be that the structure observed is the final result of a fully relaxed structure 
where high energy deformations have relaxed through formation of topological defects and the lower energy elastic 
deformations remain. Topological defects could also be preexisting in the underlying graphite material. Even if the 
observed structure in graphene is not due to these types of defects, it is worth to study how their presence would affect 
the electronic and elastic properties of the material. In terms of the elasticity theory the pentagons and heptagons 
represent disclinations of the lattice. Each disclination is made by either removing or inserting a wedge of material 
as shown in Fig. [4] for the particular case of a pentagon. Such a procedure would create an enormous elastic strain in 
fiatland, i.e., if you insist in gluing the edges in flat space. The structure is much more - or fully - relaxed if the crystal 
is allowed to bend in the third direction forming a cone (or a saddle point structure in the case of excess angle). This 
procedure is not possible in three dimensional crystals where disclinations are very rare. A dislocation can be formed 
by joining a pentagon and a heptagon through a line which will represent the Burgess vector of the dislocation. The 
Stone- Wales defects [75] play a very important role in the structural stability of the previously existing graphene 
structures. They are formed by two pent agon- heptagon pairs oriented as displayed in Fig. [Sjand are the result of a 90 
degrees rotation of a bond in the honeycomb lattice. In terms of the elasticity theory they can be represented as two 
adjacent dislocations. The spontaneous formation of the Stone- Wales defects is known to be the main mechanism to 
relief strain in plastically deformed nanotubes [78 . 

The electronic properties of topological defects in graphene have been studied at length in the literature by applying 
various numerical methods [TOHST] , The main founding is that they induce charge inhomogeneities in the samples 
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FIG. 3: A Stone- Wales defect is made by a 90 degrees rotation of a bond in the honeycomb lattice what produces two 
pentagon-heptagon rings. 




FIG. 4: (Color online) Left: Effect of a pentagonal defect in a graphene layer. Right: Cut-and-paste procedure to form the 
pentagonal defect. The points at the edges are connected by a link what induces the frustration of the bipartite character of 
the lattice at the seam (taken with permission from Ref. [88]). 



with characteristic patterns that can be observed in STM and EFM experiments. Substitution of a hexagonal ring by 
odd-sided polygons break the electron- hole symmetry which is preserved by even-sided rings. Pentagonal rings have 
been suggested as excellent field emitters [82 . In the recent field of topological insulators [56l |57] topological defects 
are able to bind excitations with fractionalized quantum numbers 23 El]- The tight-binding approach for graphene 
in the presence of pentagonal and heptagonal rings has been worked out in the context of nanotubes in the pioneer 
paper of ref. [85 . Analytical approaches give rise to the appearance of various gauge fields which we will review. An 
interesting question is whether these "plastic" deformations can be distinguished from the elastic deformations which 
also give rise to gauge fields in the continuum limit [36 . One of the main papers to address this issue that appeared 
again in the nanotube literature is [86 . It has also been discussed very recently in . 

We will first review the case of disclinations where a hexagonal ring of the honeycomb lattice is replaced by either a 
pentagon or a heptagon. We will deal with the continuum model under the assumption that the electronic structure 
of graphene is well described by (two) massless Dirac equations in two space dimensions. 

B. The gauge approach to disclinations 

We will see that modelling disclinations in the continuum gives rise to three different gauge fields coupled to the 
electronic degrees of freedom. One is associated to the phase mismatch seen by the electron when circling around 
the defect (holonomy). This one is independent of the lattice structure and will take care of the conical shape of 
the graphene electron energy surfaces [48l [89ti92] . The second class has to do with the frustration of the triangular 
sublattice ordering that occurs in the seam (see Fig. |4| and the exchange of the Fermi points that occurs when going 
around a pentagon or an heptagon. These effects where described in full detail in |48l [89] and can be modelled by 
non-Abelian gauge fields. Finally the presence of various defects forces the introduction of a third type of gauge field 
associated to the fact that the tight binding phase does not commute with the holonomy associated with the exchange 
of the Fermi points [49l [93l [94] . 
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1. Two gauge fields for a single disclination 

In relativistic quantum field theory where there is a very tight connection between the spin and the statistics, a spin 
1/2 particle is described by a field that belongs to the so-called spinor representation of the Lorentz group. A distinct 
characteristic of spinor representations is that when they move around a closed path they pick up a minus sign (upon 
a 27r rotation they acquire a phase of tt). A spinor going around a closed path encircling the pentagon in Fig. 4 will 
acquire a phase proportional to one half the total angle of the path (j) = 27r(l — 1/6). When solving the Shrodinger 
equation H"^ = with the Hamiltonian Q this condition has to be imposed on the wave function as a boundary 
condition which is often hard to deal with. A way to incorporate the constraint attached to the boundary condition 
to the Hamiltonian is to remind the Bohm-Aharonov (BA) effect and substitute the pentagon by a fictitious magnetic 
field located at the same point. The fiux of the field can be adjusted so that the phase acquired by the spinor when 
going around the vector potential is the same as the one induced by going around the pentagon. This procedure has 
the advantage of being easily generalized to conical defects of arbitrary opening angle. In the BA effect the phase is 
proportional to the circulation of the vector potential along the closed path: 

^ A.dr = 27r(l - 1/6). (25) 



The simplest vector potential having the property (25) is a vortex: 



A / N 57r / —y X \ 57r„^, , 

A{x,y) = ^ = -^V^ r , 26 

3 + + y 3 

which in polar coordinates reduces to the gradient of the polar angle 6. The presence of a pentagon in the lattice has 
an additional consequence which was discussed at length in [89 : the two sublattices of black and white points are 
exchanged by going around a conical singularity, the admissible wave functions are made of pairs of plane waves with 
opposite momenta as described in Se. In momentum space the inversion with respect to the origin exchanges the two 
independent classes of Fermi points. It becomes clear that, for the mentioned honeycomb lattices, the states of the 
theory have to accommodate into the spectrum of two coupled Dirac spinors. This new condition can be fulfilled by 
attaching a quantum number (fiavor) to the Fermi points, and considering the two bi-spinors ?/^^(r),z = +, — as the 
two components of an SU(2) fiavor doublet. The vector field will now be a non-Abelian gauge field able to rotate the 
spinors in this fiavor space. The full boundary condition to impose on the spinor when circling a pentagon (or any 
conical singularity of arbitrary defect angle (p) is 

^{0 = 0) = Tc^(<9 = 27r) ^ (27) 

^{0 = 0) = exp( i AaT''dr)^{0 = 27r), 



A.dr = 2t: — ip. 

where A^ are a set of gauge fields and a set of matrices related to the fiavor SU(2) degree of freedom of the system. 
In this way we end up with the Dirac equation coupled to a gauge potential: 

H = -ihv^j.d + ^/7^7.A(r), (28) 

where is the Fermi velocity, 7^ is a 4 x 4 matrix constructed from the Pauli matrices chosen so as to accomplish 
the gauge transformation produced by the given defect, the Latin indices run over the two spatial dimensions and g 
is a coupling parameter. The external field Ai{f) takes the form of a vortex 

^'(rl = §^^^^^^, (29) 

and the constant ^ is a parameter that represents the strength of the vortex: ^ = f Adf and is related to the opening 
angle of the defect. 

Since disclinations corresponding to even-sided polygons do not frustrate the sublattice symmetry, the associated 
matrix 7*^ in eq. (28) for this case is 



(30) 
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2. Generalization to several defects: gauge fields arising from the holonomy 

A generalization of the gauge approach to include various topological defects was presented in [49l |92l |93] . The 
strategy consists of determining the phase of the gauge field by parallel transporting the state in suitable form along 
a closed curve surrounding all the defects. 

^{0 = 0) = Tc^(<9 = 27r) (31) 
^ ^((9 = 0) = exp( l AaT''dr)^{0 = 27r), 

^ c 

where are a set of gauge fields and a set of matrices related to the pseudospin degrees of freedom of the system. 
When dealing with multiple defects, we must consider a curve surrounding all of them, as the one sketched in Fig. 

El 




FIG. 5: Scheme of a prototypical curve enclosing multiple defects in which the state will be parallel transported. 

The contour C is made of closed circles enclosing each defect and straight paths linking all the contours to a fixed 
origin. The parallel transport operator Pc associated to the closed path is thus a composition of transport operators 
over each piece: 

P = P^i-Pi- P-/ • ... • P^N • Pn ■ P~N. (32) 

where Pi is the phase associated with going around a single defect discussed in the previous subsection and P^i is the 
chiral phase acquired by the wave function while travelling in the straight paths and depends on the distances (n, m) 
between defects. 



The total holonomy turns out to be 



exp f ^ yZi^j - ^i)^3 I , (33) 




where N is the total number of defects equal to the sum of pentagons (A/'+) and heptagons (A/"-), and (n, m) is the 
position of the given defect in the usual chiral basis for the graphene lattice. Then we see that in the modelling of 
disclinations with the gauge approach described in this section there appear three different gauge fields to incorporate 
to the Dirac equation, which couple to the matrices <J3, r2 and ts: 

1. Gauge fields due to the holonomy: deal with the deficit or surplus of graphene slices 

A^ = exp(-z7rcr3/6). 

This boundary condition only affects the AB pseudo-spin structure. 

2. The second type deals with the identification of the two sublattices 

= exp(i7rr2/2). 



This part of the boundary condition affects the valley (i^, K') structure of the spinor: the two conical points are 
intertwined. 
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3. A third holonomy appears when more than one defect is considered: 

A3 = exp ^(n^- - mj)T3 

This boundary condition comes from the fact that the tight binding phase does not commute with the holonomy 
associated with r2. 



From (33) it is easy to see that for a Stone- Wehs defect we have A/'+ = N- = 2, = = 1 and = 1 and hence 
there is no gauge field associated to it. This casts some doubts on the completeness of the gauge description as it is 
known that such a defect alters the electronic structure of the graphene sheet [50l [88] . In the next section we will see 
that curvature effects account for the missing piece. 



C. Modelling curvature. Dirac fermions in curved space 



Many of the unusual properties of graphene arise from the fact that its quasiparticles are described by Dirac 
spinors. The special structure of the spinors is linked in particular to the Klein paradox and to the absence of back 
scattering in the samples [25 . When trying to study the possible consequences of the ripples found in the samples 
on the electronic properties it looks very natural to use the formalism of quantum field theory in curved spaces [95] . 
The main justification to use this approach is the robustness of the Dirac description versus deformations of the 
underlying lattice described in Section \ll B[ We must notice that, although the formalism is similar to that used in 
general relativity, here the geometry (curvature) of the space is fixed and the description is purely geometrical: the 
aim is to construct a Hamiltonian with scalar quantities made out of vectors or spinors in a curved background. This 
approach was applied to study the electronic spectrum of the fullerenes in [48l Ell [89l EH |96] . A somehow similar 
approach is used in [97] in the framework of the equivalence between the theory of defects in solids and the three- 
dimensional gravity [4^ to study the response of electromagnetic charges to conical defects in planar graphene. The 
geodesic around a dislocation has been treated in [98j. The coupling of the electronic degrees of freedom of planar 
graphene to conical defects within the geometric formalism has been explored in [50l [88l [SU [99] . There it was found 
that a distribution of pentagons and heptagons induces characteristic inhomogeneities in the density of states of the 
graphene surface. It was shown that the charge density was enhanced around the pentagonal defects and depressed 
near the heptagons, a result that was also predicted earlier in numerical calculations [79 . The effect was very localized 
and disappeared a few lattice constants away from the defect. 

In order to investigate the effect of pure curvature on the electronic properties of graphene, the conical defects 
studied previously present two difficulties. First, they correspond to surfaces with zero intrinsic curvature; moreover, 
the extrinsic curvature is accumulated at the apex of the cone where the surface has a singularity. It is then not clear 
if the results obtained - which look similar to the ones got with vacancies in [100 - are due to the singularity or to the 
curvature. To disentangle the two effects a flat graphene sheet with a smooth curved portion with intrinsic curvature 
was analyzed in |101j . In all cases the curvature gives rise to gauge flelds with a derivative coupling. In what follows 
we will sketch the formalism and give a summary of the results obtained so far within the geometric approach. A 
derivation of the main geometrical factors associated to the formalism is given in Appendix [A] 

The dynamics of a massless Dirac spinor in a curved spacetime with the metric tensor g^^ is governed by the 
modifled Dirac equation: 

i7^(r)V^V^ = 0. (34) 

The two geometric objects entering the Dirac equation are the Pauli (gamma) matrices and the derivative of the 
spinor fleld. Both are vectors that need to be properly deflned in a curved space. The curved space 7 matrices must 
obey the generalized anticommutation relations: 

{7''(r),7'^(r)} = 25'"^(r), 

and become functions dependent on the point of the space. They are related to the constant flat matrices through 
the "flelbeins" described in the Appendix |A| The covariant derivative operator is deflned as 

= - 

where F^ is the spin connection of the spinor fleld whose detailed construction is also given in the Appendix [A] 

Once the metric of the curved space is known there is a standard procedure to get the geometric factors that enter 
into the Dirac equation. In the modelling of the graphene ripples, the metric can be treated as a smooth perturbation 
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FIG. 6: (Color online) A smooth curved bump in the graphene sheet (taken with permission from Ref. [IQljl. 



of the flat surface and physical results are obtained by a kind of perturbation theory. Very often, the final result can 
be casted in the form of the flat Dirac problem in the presence of a potential induced by the curvature. In ref. |lQlj 
a smooth general metric was considered of which an example is given by the gaussian shape of Fig. [6j A simple way 
to derive the metric of a smooth shape (spherically symmetric for convenience) is as follows: Start by embedding a 
two-dimensional surface with polar symmetry in three-dimensional space (described in cylindrical coordinates). The 
surface is defined by a function z{r) giving the height with respect to the flat surface z = 0, and parametrized by the 
polar coordinates of its projection onto the z = plane. To obtain the metric we compute 

dz^ = (^^^ dr^ = af{r)dr^, (35) 

and substitute for the line element: 

ds^ = dr^ + r^dO^ + dz^ = (1 + af{r)) dr^ + r^dO^. (36) 
For the particular example of the gaussian bump given in Fig. [6j 2: = Aexp(— r^/6^), the line element reads: 

ds^ = - (1 + af) dr^ - r'^dO'^, (37) 
where a is the ratio of the height to the mean width of the gaussian that can be used as a perturbative parameter 



over the flat case. The metric (37) can be written in a more usual form: 



= ! -<^+;^M'_°. I- (38, 



Comparing the flat Hamiltonian written in polar coordinates: 



with the curved Hamiltonian 

it can be seen that the curved bump induces an effective Fermi velocity Vr in the radial direction given by 

Vr{r,e)=VF{l + af{r))-^/^ (41) 
and an effective magnetic field perpendicular to the graphene sheet given by 



14 



The magnitude of this effective magnetic field is estimated to be of the order of 0.5 to 2-3 Tesla in the region spanned 
by the bump, compatible with the estimations given in [41], and it plays the same role in the issue of the weak 
localization of graphene as the effective magnetic fields discussed there and in [42] . 

To first order in the perturbative parameter a the effective potential that appears in the Dirac equation in polar 
coordinates is 



il-dr). (43) 



The physical properties of the system are obtained from the Green's function in the curved space as explained in 
Section IVBl 

The metric approach was used to describe topological defects in ^ . The proposed metric was as a generalization 
to the one used to model cosmic strings [lQ2j : 

ds^ = -df + e-^^^^^y\dx^ + dy^), (44) 

where 

N 



A(r) = ^4/i,log(r,), 



lii is related to the defect angle of the disclination. In the case of having heptagonal defects with an "excess" angle, 
the sign of fii is negative. Here 

r, = [{x-a,f + {y-hfY'\ 

where (a^, hi) are the positions of the defects. The curved space formalism applied to this geometry induces an effective 
potential: 

V[uj, r) = 2iA7^ao + iK-f^ dj + ^7^'(5^-A), (45) 

which gives rise to spacial inhomogeneities in the local density of states as discussed in Section |V| The fictitious gauge 
field coming from the spin connection is similar to the one encountered in elasticity formalisms. A distinctive feature 
of the covariant approach is the space dependent coefficient in the kinetic term of the effective Hamiltonian that can 
be interpreted as a space-dependent Fermi velocity that could eventually be used to test the validity of the model. 
The influence of this extra term on the minimal conductivity has been analyzed in [ lQ3j . 



D. The electronic spectrum of fullerenes 



The ideas of modelling conical defects by gauge flelds and adding geometric curvature were established in the early 
publication [89 where they were applied to compute the electronic spectrum of the spherical fullerenes. The most 
popular of them, the "buckyball" Ceo has the shape of a soccer ball, consisting of 12 pentagons and 20 hexagons. The 
structure is that of a truncated icosahedron. There is a full set of fullerenes of the same symmetry consisting of adding 
crowns of hexagons around the pentagons. Following the analysis described before, modelling the Ceo molecule can 
be done with 12 non-Abelian vortices that would account for both the singular curvature, and the exchange of the 
Fermi points associated to the pentagons. Given the compact nature of the system, the model proposed in |89j was 
to solve the Dirac equation on the surface of a sphere with a flctitious magnetic monopole at its center. The sphere 
smears the curvature accumulated at each pentagon, and the charge g of the flctitious magnetic monopole is adjusted 
by adding up the individual fluxes of all the lines: 

i=l 

N being the number of conical singularities on the surface. It is interesting to note that the value of g required for 
the icosahedron {N = 12), g = 3/2, is compatible with the standard quantization condition of the monopole charge. 

The given model reproduces the low-energy spectra and correct numbers of zero modes for the Ceo family of 
fullerenes and works better in the limit of large number of points. The spectrum is obtained by solving the eigenvalue 
problem for the covariant Dirac operator 



icr'^e^(V^ - iA^)^n = en^n^ a, /i = 1, 2, 



(47) 
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where ej^ is the zweibein for the sphere, and the geometrical factors in spherical coordinates are 



V(/> = 0(1) -[(J-^cr^l COS I* 



= 0, 



^0 = g cos Or^'^\ 



The Dirac equation becomes 



R 
R 



1 



sin(l9) 
1 



2cos(i9) 



sin(l9) ^ 2cos(l9) 



2sin(l9) 
2sin(l9) 



(48) 



where R is the radius of the sphere. 

The Dirac equation can be diagonalized by introducing a generalized angular momentum operator J taking into 
account spinor indices, magnetic field and curvature. We find 



J- = -e 

Jz 



1 



sin(l9) 



2sin((9) 



+ 



sin((9) ' 



> 9-1 
2sin(l9) 



-ids 



By squaring the Dirac operator each of the spinor components obeys the equation 

2 



(49) 



(50) 



where J is the eigenvalue of the angular momentum, which is an integer. The degeneracy of a state is 2J + 1. From 
Eq. (50) it can be seen that there is minimum value of the angular momentum J dictated by ^, so that the number 
of zero modes in the spectrum depends exclusively on the value of the monopole charge, i. e. on the number of 
pentagons in the lattice. It is also interesting to note that, since the sphere has a constant curvature, there are no 
inhomogeneous geometrical factors and all the effect of the curvature reduces to a rescaling of the energy. Substituting 
g = 3/2 in Eq. (50) gives the spectrum of the Ceo family of fullerenes: the ground state is a couple of triplets lying 
at zero energy, and the first excited levels follow the degeneracy of the sphere up to the point in which the highest 
dimension of an irreducible representation of the icosahedron symmetry group is reached. The model works very good 
given its simplicity. The elastic properties of the fullerenes were analyzed within this approach in [48 . 



E. Dislocations and torsion: a general relativity approach 



The gravity connection can be pushed forward and used to model a density of dislocations in the space in the 
continuum limit by adding torsion to the space spanned by the lattice. The connection of torsion with the continuum 
theory of crystal dislocations goes back to Kondo [104 and has been formalized in (SHIl^. Various aspects of the 
problem have been explored in ^105j . A nice review on the relation of gravity with topological defects in solids is 
|106j . The geometric approach to defects in solids [36l I107j relates the metric of the curved crystalline surface with 
the deformation tensor and establishes that disclinations are associated to the curvature tensor and a finite density 
of dislocations generates a torsion term. 

Dislocations in graphene are made of pentagon-heptagon pairs and they have been widely studied in connection 
with the properties of carbon nanotubes [74 and, more recently, in the fiat graphene sheets or ribbons fl08 | 1109] , 
The covariant approach described above was extended in |110j to include a connection with torsion. As we discussed 
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in Section II D the minimal coupling of any geometrical or physical fields to the Dirac spinors adopts always the form 



of a covariant derivative: 

d^^D^ = d^ + A^, (51) 

where the given vector can be an electromagnetic potential induced by a real electromagnetic field or any other 
real or fictitious gauge field associated to deformations or to geometrical factors. In the case of having a density of 
dislocations in the graphene sheet modelled by torsion we can construct two potential candidates to gauge fields: 

~ 9 ^^upfi -) Sfi — ^ pu pa^u pa •) (^^) 

where T^^p is the rank three torsion tensor related to the the antisymmetric part of the connection: 

The field is an ordinary vector and can be associated to the density of edge dislocations while is an axial 
(pseudo) vector associated to the density of screw dislocations. 

IV. STRAIN FIELDS 
A. Modulations of the tight binding parameters 

In this section we will consider gauge fields which arise in defect-free graphene due to smooth elastic deformations. 
The latter can be created by applying external stress to graphene flakes. It is important to notice that these defor- 
mations and the related gauge fields are also intrinsic for graphene at finite temperatures due to thermal fluctuations 



(see Section IV D). These gauge fields were discussed in a context of electron-phonon interactions in nanotubes prior 
to the synthesis of graphene in [62l llllj . 

Let us start with a tight-binding model of the electronic structure of graphene in the nearest-neighbor approximation. 
In the deformed graphene all bonds are, in general, nonequivalent and the three nearest-neighbors hopping parameters 
ti can be all different. If one repeats the derivation of the Dirac Hamiltonian in the effective mass approximation 
but with nonequal hopping parameters one finds for the electron states in the vicinity of the K point the effective 
Hamiltonian: 

H = -ihvra (v - il^ , (54) 

where [13 |62l [m] 

/3 

A, = ^(^3-^2), 

Ay = ^(t2 + t3-2ti), (55) 

(note that, following Ref. [6^ and most of the subsequent works, we use a choice of the axes x, y different from that 
made in Refs. [271 111!")- In the weakly deformed lattice, assuming that the atomic displacements u are small in 
comparison with the lattice constant a one has 

ti=t^ ^Pi {ui - uo) , (56) 
where pi are the nearest-neighbor vectors, uq is the displacement vector for the central atom, and 

/^ = -|j^-2 (57) 
am a 

is the electron Griineisen parameter (for more details, see below). The continuum limit (elasticity theory) with a 
displacement field u (r) is performed by making the substitution 

Ui-uo (X (pi V) u (r) , (58) 
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from where we obtain the effective gauge field [HH ESJ: 

— ^ CI ^yy ) ' 

Ay —C—^Uxy^ (^9) 

where c is a numerical factor depending on the detailed model of chemical bonding; in what follows we will make our 
estimations with the value c = 1. 



For the other valley, K' ^ the sign of the vector potential (59) is opposite, in agreement with the requirement of 
time-reversal invariance. Indeed, atomic displacements can break effectively the time-reversal symmetry (involving 
inversion of the wave vector) for a given valley but not for the crystal in general. 

If the deformation is not a pure shear and the dilatation does not obey the condition Vu = Uxx + Uyy = 0, a scalar 
potential proportional to Vu also arises [62]. The gauge field ([59| is proportional to the deformation tensor which 
is directly involved in the density of elastic energy (see the next Section). This means that the problem as a whole 
is not gauge invariant. For example we could add to the gauge potential the gradient of a scalar function and shift 



the electronic wave function by a phase as discussed in Section II C leaving the electronic part of the Lagrangian 



invariant but these transformations will change the elastic part. The simplest argument to understand this lack of 



gauge invariance is to notice that the kinetic energy term for the elastic part is isotropic (see Eq. (62)) while the 



effective vector field given by Eq. (59) is not. 



On general symmetry grounds, strains can also lead to scalar potential [62l |63], 

V(r) = g{uxx + Uyy) (60) 

This term enters in the diagonal elements of the Dirac equation. The most recent calculations suggest that g « 4eV for 
a neutral graphene single layer [1 12 , although higher values, g ^ 20eV have been cited in earlier literature [62l 11131 1114j . 
The parameters j3 and g describe the coupling to acoustical in graphene, and they are usually considered simultane- 
ously, as they both contribute to the deformation potential, D. Combinations of the type = ci x {v'^^'^) / o? x 
where Ci and C2 are numerical constants determine the phonon contribution to transport coefficients, like the resistiv- 
ity. Microscopic models for the deformation potential have considered only the effect of the gauge potential |115l4118j . 
or both the scalar and gauge potentials |119j . 



B. Gauge fields as function of the in plane strains 

Long wavelength distortions within the graphene layers can be described using the strain tensor [120], Uo,^. This 
tensor defines the local deformation of the lattice. Given the deformation field u{x^ y)^ the strain tensor in the linear 
approximation is: 



Ua^ = . (61) 



The elastic energy is: 



Eeias = \ j d^^i^xx + Uyyf ^ (1 j d?f[ul^ + v^y + 2uly) , (62) 



where A and ji are the elastic Lame coefficients. We can also define the stress tensor, cFcxj3 = dEeias/dua(3^ 

^xx — 

A {UxX + Uyy) + 2llU 

XX 

^yy — ^ {Uxx H~ Uyy) -\- IflUyy 

(Jxy = 2llUxy. (63) 



As discussed in Section II D[ strains modify the bond lengths and the hopping between tt orbitals. The dependence of 



the hopping on the interatomic distance is approximately the same in the organic compounds with sp^ coordination. 



and it is parametrized by the dimensionless constant /3 (57). A change in the hopping between nearest neighbor 



atoms modifies the terms in the Hamiltonian which couple the two sublattices. In the previous section this was done 



explicitly, in the nearest-neighbor approximation. Actually, the expression (59) is more general. In the continuum 
limit, this coupling must take the form of a gauge field with the matrix structure of the Pauli matrices, ai and (72. 
A vector can be constructed from the strain tensor by contracting it with a rank three tensor. This tensor should be 
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invariant under the discrete set of symmetries of the honeycomb lattice. The only possible tensor with these properties 
has the non zero components: 

K — 1 

Kyyx Kyxy ^xyy 1? (^'^) 

where the x axis coincides with one of the unit vectors of the lattice. The same tensor describes the first deviation 
from complete isotropy in the bands of the electrons in the honeycomb lattice, the so called trigonal warping [19 . 



From Eq. (64) t can be seen that the gauge field must take the form (|59|). 



The combination of strains which enter in the gauge field, Eq. (59), can also be written in terms of the stress tensor. 



2/i 
2/i ' 



2u,y = (65) 



where /i is a Lame coefficient (shear modulus). 

The calculation of the strains in a two dimensional system with only in-plane displacements is drastically simplified 
in comparison with a generic three dimensional case, as the stresses can be written in terms of a function [120 : 

/(x, y) = Re [{x - iy)g{x + iy) + h{x + iy)] , (66) 

where g{z) and h{z) are analytic functions of z = x -\- iy. The stress tensor is: 

_ 

- Qy2 

""yy ~ dx^ 

A pure shear deformation has = ~^yyi which implies that g{z) = 0. Then we get the effective gauge potential: 

ocRe [h"{z)] 

Ay 0Qlm[h"{z)], (68) 

and the effective magnetic field is: 

B = -dyAx + dxAy (X Im [h''\z)] . (69) 

A pure shear deformation that leads to a constant effective magnetic field is given by h{z) = Az^^ where A is a 
constant. A general deformation which satisfies the equilibrium equations of elasticity and which leads to a constant 
effective gauge field is determined by the function f{z) = Az^ + Bzz^. 

A constant effective magnetic field requires a strain tensor which increases linearly with position, Uij = u{fn)/L, 
where n is a given unit vector, and L is the size of the system. The maximum stresses are of order u. The effective 
magnetic field has an associated magnetic length: 

where a is the lattice constant. For u = 0.01 and L = 1/im we obtain Ib ~ lO^nm, which corresponds to a magnetic 
field of 0.5T. A description of the strains and effective fields induced by displacements with radial symmetry is given 
in the Appendix [B| 



C. Out of plane displacements 



When the layer can be distorted along the out of plane direction, the strain tensor becomes: 
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FIG. 7: Sketch of a suspended rectangular flake. 



where h is the displacement along the third dimension. The distortion of the bonds is proportional to the square of h. 

The simplest and most ubiquitous deformation expected in suspended graphene is due to the force induced by the 
electric field between the gate and the flake [121 . The total energy of a monolayer flake under a constant vertical 
force, as function of the maximum vertical deformation, /i, can be written as |12Qj : 

Etot = Ebend + Eelas + E field ^ C^end^-^ + {^xX + C^/i) j-^ + CfieldJ^hL , (72) 

where Cbend^cx^Cfi and Cfieid are numerical constants, k, ^ 1 eY is the bending rigidity of the flake, A ^ 2.4 eV 
and /i ^ 9.9 eV are the in plane elastic Lame coefficients and = 2Tie^in? is applied field, and n is the 
carrier density that it induces. The bending energy in Eq. ( 72 ) is negligible with respect to the in plane elastic energy 
for h > ^/hi/m8ix{X^ ji) ^ 1 — 3 A. Neglecting it, we obtain: 

h^CLi A , (73) 



max(A, /i) 

where C is a numerical constant which depends on details of the shape of the flake, and the ratio X/ ja. The neglect 
of the bending energy is justified for n^L"^ > y^^^Vl^^^^'^^^l^T/^ ~ 0.1. As = nL^ is the total number of carriers in 
the flake, the bending energy can be neglected in all realistic situations. For n ~ lO^^cm"^ and L ~ 1/im, we obtain 
h ~ 30nm. 



The strains associated to the deformation in Eq. ( 73 ) are of order: 



2/3 



max(A, /i) 

leading to an effective magnetic field characterized by a magnetic length, Ib- 



^ {Vhf ^ C ( ) , (74) 



1 C / e^n^L 



2/3 

ii - /?""^ [:^^)) ' 

where a is the lattice constant. For n ~ lO^^cm"^ and L ^ 1/im, we obtain ~ 30 — lOOnm, which corresponds to 
effective fields of 2 - lOT. 

The effective magnetic field is determined by the full strain tensor, which includes a contribution from the in plane 
displacements |122j . If the height profile, /i(r), is known, we can define: 

dh dh 

Then, the Fourier component of the effective field, including the effect of the in plane strains can be written as: 

B(k) (X iky — ^ kyfxx{k) + k'^fyy{k) — 2kxkyfxy{k) , (77) 

(A + 2/ij|A:|^ >- 

where k is the wave vector. These corrections do not change the earlier order of magnitude estimates for the effective 
field, except in cases with high symmetry. For instance, when h is constant along one direction, i. e. the relaxation 
of the in plane displacements lead to a stress tensor such that ctxx = constant, (jyy = dxy = 0, and to a constant 
effective gauge field. 
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FIG. 8: Transmission through a suspended rectangular flake as function of incident angle. The carrier density is p = 10 cm~ , 
the length of the flake is L = lOOnm, and the maximum deformation is Hq — 3nm. Red and green lines correspond to the two 
valleys in the Brillouin Zone (see text for details). 



FIG. 9: Sketch of a flexible membrane (solid line), h is the out of plane deviation with respect to the z = plane (dashed line) 
defined by the center of mass. The unit vector n and no are the normals to each point in the membrane and in the reference 
plane respectively. 



From the general analysis of the previous sections it follows that corrugations of graphene with spatial scale much 
larger than the interatomic distance a lead to the formation of an Abelian gauge field. It turns out that for two- 
dimensional crystals at finite temperatures the corrugations are unavoidable, due to the instability of the crystal with 
respect to the bending fluctuations. In this section we will review the corresponding theoretical results. 

The standard theory of lattice dynamics and thermodynamics |123j starts with the assumption that typical ampli- 
tudes of atomic displacements u due to thermal fluctuations are small in comparison with the interatomic distances, 
(v?) <C a^. It was understood as early as in 1930th p!| I124fll26] that for two-dimensional crystals at finite tempera- 
tures this cannot be the case, due to the logarithmic divergence of {u^)- This has given rise to the hypothesis that 
long-range crystal order cannot exist in two dimensions which has been proven later, as a particular case of Mermin- 
Wagner theorem |127j . Therefore the experimental realization of graphene and other truly two-dimensional crystals 
p!2l [T3] was really surprising. The point is that these objects are not two-dimensional crystals in two-dimensional 
space but two-dimensional crystals in three-dimensional space; it follows from a general theory of flexible membranes 
[128 that such systems can exist but cannot be really flat [T7l[43]. Very soon after the synthesis of graphene corru- 
gations of freely suspended graphene membrane were observed experimentally [43j; the existence of intrinsic ripples 
in graphene due to thermal instability has also been confirmed by atomistic Monte Carlo simulations [129 . 

Let us start with the model of a continuum infinitely thin membrane (Fig. |9|. It is characterized by the dependence 
of out-of-plane deformation on coordinates, h{x^y). The normal vector to the membrane at a given point has 
components [128j 




D. Thermal fluctuations and intrinsic ripples in graphene 




(78) 



For liquid membranes, according to the Helfrich model [130\ I131j the density of elastic energy contains two terms. 
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FIG. 10: A representative atomic configuration in Monte Carlo simulations 
(courtesy A.Fasolino). 



[T29l [T33] for graphene at room temperature 



proportional to the square of mean curvature H'^ and to the Gaussian curvature K but the second one gives just a 
constant contribution to the total energy (due to Gauss-Bonnet theorem) and is irrelevant for our analysis. In the 
lowest-order (harmonic) approximation H (xV^h and, thus, the bending energy has the form 



(79) 



In the harmonic approximation, the correlation function of the Fourier components of /i(r), \ \hc^\ j is calculated 
immediately: 



A ' 



Thus, the mean square amplitude of out-of-plane displacement is 

,2\ T. 



(80) 



(81) 



where L is the sample size and we have introduced a cutoff of divergent integral in Eq. (81) at 



- l/L. Thus, 

typical amplitude of out-of-plane fluctuations is proportional to the sample size. The correlation function of the 
Fourier components of the normal vectors. 



(82) 



is also singular at small wave vectors and, thus, the correlator (n(^R^ ^(0)^ behaves logarithmically at ^ oo 

whereas for a flat membrane it should tend to a constant (normals are more or less parallel at large distances). As a 
result, the model with the Hamiltonian (79) describes not a flat but a crumpled membrane. 



The higher-order terms in the expansion of the mean curvature H m\/h renormalize the effective bending rigidity 

[Ml, 



3T L 

«e//(T)=«--ln-, 



(83) 



thus, the membrane becomes less rigid as the temperature increases. 

The situation is totally different for crystalline membranes where the coupling between out-of-plane (bending) h (r) 
and in-plane u (r) atomic displacement becomes crucially important [128 , 134 . It follows from the elasticity theory 
|12Q[ I135j that the total elastic energy in this case is written as (cf. Eq. ([62|) 



n 



A 



2 ^OiOi 



(84) 



Uoip is given by Eq. (71) We keep there quadratic terms in the derivatives of out-of-plane deformation but not in the 



derivatives of in-plane deformations since the former are, at the average, much larger. The model described by the 
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FIG. 11: Portion of one typical atomic configuration of graphene sample at room temperature. The numbers indicate the bond 
length in A (taken with permission from Ref. [129 ). 



Hamiltonian (84) is essentially nonlinear. Since the Hamiltonian is quadratic in in-plane components of the distortion 



tensor, the partition function can be integrated over the latter and the effective Hamiltonian for out-of-plane 

ox jj 

deformations only can be exactly derived |128[ 11341 11361 1137j : 



n 



I 



Ko 



dh dh 



(85) 



where Pj^ = 5^/3 — {d'^ / dxo^dx^^ /V^ is the operator separating transverse components and 



4.11 [ii + A) 
2/i + A 



(86) 



The Hamiltonian (85) describes interactions of soft (long- wavelength) fluctuations, similar to the Ginzburg-Landau- 



Wilson Hamiltonian for the critical point [6l |7l [138 . The difference is that two-dimensional systems are "critical" 
at any finite temperatures. In analogy with the theory of critical phenomena, one can introduce |128[ I134j scaling 
hypothesis for correlation functions and the corresponding exponents 77,7?^C related with the behavior of effective 
bending rigidity, Keff {q) oc elastic moduli, jHeff (q) , Xeff {q) q^ and characteristic out-of-plane deformation, 



(87) 



{h)'^ oc L^. All these exponent can be expressed through the only exponent r]: 

C = l-7?/2, r,' = 2{l-ri) 
The normal-normal correlation function, instead of Eq.(l82l), behaves as 



G{q) (X q 



-(2-r7) 



at g ^ 0, and the membrane is "flat" (in a sense that (n (^R^ ^(0)^ ^ const at R ^ oo) for < < 1. However, 
the amplitude of out-of-plane displacements grows with the sample size which means that the membrane is essentially 

corrugated, in contrast with the surface of three-dimensional crystals where \J {h)'^ remains much smaller than the 
interatomic distances. Note that the in-plane phonon modes become even softer due to the coupling with the out- 
of-plane ones, since /ig// {q) i^eff {q) ^ ^ ai q ^ 0. In this sense, there is no real crystalhne order since there is 
no infinitely narrow delta-functional Bragg peaks but, instead, some power-law singularities at the reciprocal lattice 
points so the Mermin theorem is, of course, not violated. This broadening of the Bragg peaks is clearly seen both in 
experiments [43 and simulations [47] for graphene. The continuum model (84), which is called the model of phantom 
membranes, has a transition to a crumpled phase at a temperature of the order of The term "phantom" means 
that the model does not include self- avoidance, the natural condition of true physical systems. It is assumed that 
self-avoidance removes the phase transition to the high temperature crumpled phase while the scaling properties of 
the "flat" phase remain the same as in phantom membranes |128j . To calculate the exponent 77 some approximations 
should be done. The e-expansion in e = 4 — where d is the space dimensionality, so efficient in the theory of 
critical phenomena [H [3 I138j is practically useless for low-dimensional systems. Instead, one can use an expansion 
in inverse number of component n of the vector of out-of-plane displacement h. Of course, in real physical problem 
n = 1. However, if one consider D-dimensional membrane in (i-dimensional space, n = d — D can be formally 
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FIG. 12: Radial distribution function at T =300 K and T =3,500 K as a function of interatomic distance. The arrows indicate 
the length of double (r =1.31 A), conjugated (r =1.42 A) and single (r =1.54 A) bonds (taken with permission from Ref. [129] ). 



considered as a large parameter in the problem [134, than, rj ~ 1/n. More subtle use of this parameter is the Self 
Consistent Screening Approximation |136[ 1137] which neglects vertex corrections when consider diagrammatically the 
interaction of fluctuations; the vertex corrections have formal smallness in 1/n. This approximation (for the real case 
d = 3, D = 2) yields r] = 2 (>/T5 — l) /7 ~ 0.821. A very recent non-perturbative renormalization group approach 
|139l yields a close value rj ^ 0.849. The discretized version of the model ([84| was investigated by Bowick et al |138j 
by means of Monte Carlo simulations giving r] ~ 0.72. 

For the case of crystalline membranes, contrary to the case of liquid ones, the effective bending rigidity grows with 
the temperature [128 , due to anharmonic coupling of out-of-plane and in-plane displacements. This result has a very 
simple physical meaning: it is known from the elasticity theory that corrugations strengthen plates making them 
effectively thicker |120[ 1135] , and amplitudes of corrugations increase with the temperature increase. 

The long-range orientational ( "hexatic" ) or translational (crystalline) order can be in principle destroyed by sponta- 
neous creation of topological defects, disclinations and dislocations, respectively. Whereas for two-dimensional crystals 
in two-dimensional space the disclination energy grows with the sample size as (and the dislocation one as InL) 
for the case under consideration (two-dimensional crystals in three-dimensional space) the elastic energy strongly de- 
creases due to screening by bending deformations; as a result, the disclination energy grows as InL and the dislocation 
energy remains constant in the thermodynamic limit L ^ oo |128j . This means that, strictly speaking, the long-range 
crystalline order is destroyed also via the creation of dislocations by thermal fluctuations. However, for graphene with 
its extremely strong chemical bonding (e.g., vacancy formation energy about 7 eV jl40j ) the concentration of such 
thermal dislocations should be very small; no evidences of their formation were found in atomistic simulations for 
temperatures up to 3000 K [66lfT29] . 

A quantitative information about intrinsic corrugations in graphene can be obtained using atomistic simulations 
with some interatomic interaction potential for carbon (of course, it cannot be just pairwise potential, it should be 
dependent on the angles between bonds which is crucially important for a covalent chemical bonding). The Monte 
Carlo simulations j66lll29[ri33j used the so called bond-ordered potential LCBOPII |141j . This potential is based on a 
large database of experimental and theoretical data for molecules and solids and has been proved to describe very well 
thermodynamic and structural properties of all phases of carbon and its phase diagram in a wide range of temperatures 
and pressures |14H I142j . Of particular importance here, is that bond order potentials correlate coordination to bond 
strength, allowing changes between single, double and conjugated bonds with the correct energetics. 

A typical configuration of graphene at room temperature obtained in the simulations is shown in Figure [lO] For the 
crystallite of about 10000 atoms a typical size of the height fluctuations was 0.07 nm (a half of interatomic distance) 
with a typical spatial scale 5 to 10 nm. As a result, the bond lengths are rather broadly distributed, changing between 
lengths of double and single carbon-carbon bonds (Figs. 
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and 12L The normal-normal correlation function shown 



2 from [133 ) demonstrates a crossover from the harmonic behavior (82) to anharmonic one, (88), 



[131 [139]. 



2 nm , with rj ^ 0.85. This value is in a very good agreement with the predictions of the continuum model 
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FIG. 13: Normal-normal correlation functions calculated for three systems with N= 12096, 19504, and N=37888 carbon atoms. 
The dashed lines show the asymptotic harmonic behavior with power laws with between these limiting cases (taken with 
permission from Ref. [133 J. 



V. OBSERVABLE CONSEQUENCES OF THE GAUGE FIELDS: MICROSCOPIC EFFECTS 



A. Index theorems and zero energy states 



The energy spectrum of massless Dirac fermions coupled to gauge fields has a very special property, the existence 
of zero-energy chiral modes equally shared by electrons and holes. This is a consequence of one of the most important 
theorems of modern mathematics, the Atiyah-Singer index theorem [143, 144 . This theorem has important appli- 
cations in quantum field and superstring theories |145[ I146j . It was mentioned already in one of the first works on 
graphene [14J (see also Ref. [17 ) that the anomalous ("half-integer") quantum Hall effect in single-layer graphene is 
a consequence of the Atiyah-Singer theorem (later jl47j a similar statement has been also proven for the anomalous 
quantum Hall effect in graphene bilayer [148 ). The existence of the zero-energy states for massless Dirac fermions 
in inhomogeneous magnetic fields in two dimensions was demonstrated explicitly long ago by Aharonov and Casher 
|149 L Here we present their result and discuss its relevance for graphene physics. 

Let A{x^y) is an arbitrary two-dimensional vector potential satisfying the condition VA = 0. Thus, one can 
introduce a scalar "potential" (j) (x, y) such that, that 

A^ = -dy(^,Ay=d^^ (89) 

and 

V^(^ = B (90) 

is the magnetic field. For zero-energy the Dirac equation is split into two independent equations for the spinor 
components ?/^cr (c" = ±) 

(S^ + icjdy - iA:, + a Ay) = (91) 
The potential (j) can be excluded by the substitution 

= e-^^U (92) 

where 

{d:,^iady)U=0 (93) 



and, thus, /+ and /_ are analytic and complex conjugated analytic entire functions of z = x -\- iy^ respectively. 
Since the Green's function of the Laplace operator in two dimensions 
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FIG. 14: Low energy states induced by a one dimensional ripple modelled by a sine modulation of the atomic displacements. 
The average hopping is ty = 3eV. The width of the ripple is 1200a = 168nm. The modulations of the hoppings are: Left, 
8t/t — 0, center, 5t/t — 0.02, right, 5t/t — 0.04 (taken with permission from Ref. ^152^ ). 



the solution of Eq.(90) at large distances has the asymptotics), with 



^W^^m^ (95) 



where ^ = / dxdyB is the total magnetic flux through the system. Thus, at large 

cr<3E>/27r 



(96) 



Since the entire function cannot go to zero in all directions at infinity, ipa- can be normalizable only assuming that 
cr^ > 0, that is, zero-energy solutions can exist only for one (pseudo)spin directions, depending on the sign of the 
total magnetic flux. 



Let us count how many independent solutions of Eq. (93) we have. As a basis, we can choose just polynoms 
searching the solutions of the form 

^+ = z^'e-^ (97) 



(to be specific, we consider the case $ > 0, where j = 0, 1,2... One can easily see from Eq.(96) that the solution is 
integrable with the square only assuming that j < where N is the integer part of ^/27r. 

According to the Atiyah-Singer theorem, the index of Dirac operator in our situation, that is, the difference between 
numbers of solutions with zero energy and positive and negative sigmas is equal to N. The Aharonov-Casher procedure 
gives us an explicit shape of these solutions and proves that, actually, there are solutions only for one spin projection. 
These solutions are robust (topologically protected) in a sense that they are not shifted from zero energy of broadened 
by any inhomogeneities of the magnetic field. 

It has been proven in Ref. [147 that the Atiyah-Singer theorem can be applied also for the case of bilayer graphene; 
in this situation the number of zero-energy solutions is twice the number of the solutions for single layer with the 
same magnetic flux. Explicitly, these solutions have been constructed, as a generalization of the Aharonov-Casher 
solutions, in Ref. jl50j. Thus, the zero-energy states in the case of bilayer are also topologically protected and cannot 
be broadened by any inhomogeneities of the magnetic field. 

This conclusion leads to an important consequence for the quantum Hall effect in single- and bilayer graphene. 
Whereas smooth inhomogeneities of electrostatic potential lead to a broadening of all Landau levels, like in conven- 
tional homogeneous two dimensional electron gas, and zero-energy Landau level is not an exception random 
vector potential due to ripples cannot broaden the zero-energy Landau level which, thus, should be essentially nar- 
rower than the other Landau levels. This conclusion seems to be in an agreement with experimental data on quantum 
Hall activation gaps in single-layer graphene [151 . It would be very interesting to check it also for the case of bilayer. 
If the pseudomagnetic field created by ripples is strong enough it can result in the appearance of the mid-gap states 



Figure 14 152 shows energy spectrum of a graphene ribbon with the simplest one-dimensional sine modulation 
of the atomic displacements and, thus, vector potential. One can clearly see an appearance of the midgap states 
for strong enough modulation. If the modulated scalar potential is also present the Atiyah-Singer theorem is not 
applicable and, in general, a gap can open. General conditions for opening the gap in oscillating electric and magnetic 
fields are discussed in Ref. |153j . The existence of the midgap states for sine modulation has been confirmed by ab 



initio calculations |154j . 

The midgap states are probably essential in the chemistry of graphene. It was demonstrated by a straightforward 
ab initio calculation [155 that the energy of chemisorption of hydrogen on rippled graphene is much lower than that 
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FIG. 15: (Color online)Left: Effect of the curved bump of fig. |6]on the local density of states of the graphene sheet. The color 
code is indicated in the figure. Darker (lighter) areas represent negative (positive) corrections to the density of states of the 
flat graphene sheet. Right: Correction to the density of states (dotted line) in arbitrary units, versus the shape of defect (red 
- upper - line) and curvature of the defect (yellow -middle- line) for a gaussian bump of an average width of 50 A (taken with 
permission from Ref. [101 ). 



on the flat samples indicating that the ripples can be centers of chemical activity. This effect is especially strong if 
the corrugations are large enough to create the midgap states: hydrogen leads to the splitting of the peak at zero 
energy and thus to an essential total energy decrease. This splitting results from a resonance between the midgap 
states created by hydrogen itself and those created by the ripples. 



B. Inhomogeneities in the electronic density 

One of the effects of having samples locally curved or rippled is to induce charge inhomogeneities in the samples. 
This effects were studied within the geometric (covariant) approach in Refs. [50l llOlj by computing the local density 
of states through the Green's function. In a curved surface, the spinor Green's function obeys the equation 

irej" {d^ - n^) G(r,r') = S{r - r'){-g)-i, (98) 

where the geometric factors are described in Appendix [Aj From this equation the local density of states can be 
obtained through the expression: 

p{E,r) = --ImTr[G{E, r,r)j\ (99) 

TT 

The gaussian shape of Fig. [6] was shown in |101j to induce oscillations in the density of states of the type depicted in 
Fig. [15] a). It is interesting to note that the local maximum of correction to the flat DOS corresponds to the regions 



in space where the curvature changes sign as it can be seen in Fig. 15 b). An interesting prediction of this model is 



a space variation of the Fermi velocity Vr in the radial direction given by 



vr{r.0) = ^^^= (100) 



where /(r) is the equation of the surface (35 ). It is worth mentioning that irrespective of the shape of the corrugations, 
the effective Fermi velocity in the presence of intrinsic curvature will always be less than the flat value vf- 

The same formalism applied to a given distribution of topological defects gives rise also to inhomogeneities of the 
local DOS where a local charge density is accumulated around the pentagons and repelled near the heptagons [50] . 



C. Strains, gauge fields, and weak (ant i) localization 

In normal metals, the wavevector of the electrons at the Fermi energy is comparable to the lattice spacing. This 
allows us to define wavepackets which are much smaller than the distance over which the electrons move between 
collisions, the mean free path. The trajectories of these wavepackets in an applied electric field are described by the 
classical equations of motion. 

A full analysis of the conductance of a metal requires the analysis of the effects associated to the wave nature 
of the electrons jl56| I157j . These corrections to the classical description can be studied as a series in powers of 
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FIG. 16: (Color online) Counterpropagating paths which contribute to the magnetoresistance in weak magnetic fields, a) Paths 
where the wavefunctions belong to the same valley. Although these paths are not related by time reversal symmetry, near the 
Dirac point are sufficiently similar to lead to interference effects, b) Paths related by time reversal symmetry. Interference 
between such paths is the leading contribution to the magnetoresistance if intervalley scattering is sufficiently strong. 



kpi^ where /c^ is the inverse Fermi wavelength, and i is the mean free path. The classical description is obtained 
in the limit kpi oo. The leading quantum effects arise from the interference between two counterpropagating 



paths, as schematically depicted in Fig. 16 In a normal metal, such paths show a constructive interference. As a 
result, the probability that the electron is scattered backwards is enhanced with respect to the classical case. This 
interference is suppressed in the presence of a magnetic field, and this leads to an increase of the conductivity in weak 
magnetic fields, the so called weak localization effect. When kpi 1, interference effects become strong, and they 
lead to the formation of localized wavefunctions 1 1581 1159] . and to insulating behavior. In systems with large spin-orbit 
coupling, the spin is pinned to the momentum the electron. The rotation of the spin when the electron moves around 
a closed loop adds a phase of tt to the wavefunction. As a result, closed loops show negative interference, and the 
magnetoresistance is negative, the so called weak antilocalization. 

Electrons in graphene have the sublattice and the valley degrees of freedom, which can play a role similar to the spin 
in systems with a strong spin-orbit coupling. The situation is schematically shown in FigJTBj Counterpropagating 
paths can exist where the electron resides always in the same valley, as shown in Fig. |16| a). The two paths are 
not exactly the same, because they are not related by time reversal symmetry, which exchanges the valleys. The 
equivalence between the two paths arising from the same valley becomes exact at the Dirac point. As a result, these 
paths can contribute to the magnetoresistance for weak magnetic fields, provided that intervalley scattering can be 
neglected. The contribution of these paths leads to weak antilocalization, because of the contribution to the total 
phase of the sublattice polarization, the pseudospin. This situation can be also described as the result of the existence 
of an approximate symmetry which plays the role of time reversal symmetry within each valley. 

The contribution of the paths in Fig. [16] a) is suppressed by the gauge fields due to lattice deformations [4T1|42]. The 
interference between these paths is also reduced by the inequivalence of the wavefunctions from the same valley with 
opposite momenta, which increases away from the Dirac energy, the so called trigonal warping [160 . As a result, there 
is a regime in graphene, when intervalley scattering is small, but gauge fields are significant, where there are neither 
weak localization nor weak antilocalization effects. This regime describes well the experimental situation [41" 161]. 
The absence of any magnetoresistance in the most recent experiments [161 is remarkable, because the intravalley 
scattering length associated with gauge fields or other perturbations which break the effective time reversal symmetry 
at each valley is much shorter than other scales, and it is difficult to explain by the height of the ripples observed by 
AFM measurements. 

In the presence of strong intravalley scattering, the leading paths which contribute to the quantum corrections 
of the conductivity are shown in Fig. 16 b). These are the same processes which describe the magnetoresistance in 
ordinary metals, and they lead to weak localization, and a negative magnetoresistance. Bilayer graphene, were smaller 
corrugations are expected, shows weak localization effects at low fields [162] . 
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D. Modelling disorder in graphene with random gauge fields 

The influence of ripples and other types of disorder on the transport properties of a system requires to assume a 
density of defects with some statistical distribution and to average over defects. The standard techniques of disordered 
electrons ^163j can be applied to rippled graphene by averaging over the random effective gauge fields induced by 
curvature or elastic deformations. Modelling disorder with random magnetic fields was first proposed in the context 
of disordered degenerate semiconductors ^164| 1165] . Dirac fermions in random magnetic fields were used to study 
critical point transitions between integer and fractional fillings in the Quantum Hall Effect |166[ I167j . Most of these 
works involved disordered systems without interactions what reduces effectively the dimensionality of the problem 
to two dimensions allowing the use of very powerful tools as bosonization or conformal field theory techniques [168J. 
A complete classification of disordered graphene without interaction in the random matrices language was worked 
out in |169H171] . Combined disorder and interactions was first addressed in the context of critical points between 
integer and fractional fillings in the Quantum Hall Effect |166[ \167\ I172j and later applied to the graphene problem in 
[173-175 . A recent complete survey can be found in [iTQ and the included references. A computation of the effect of 
the curvature induced by topological defects on the conductivity of neutral graphene using the random field approach 
was done in [103 . 

The low energy excitations in graphene combining the two Dirac points can be described by the four-dimensional 
Hamiltonian: 

Ho = ivF J dPx^{x)^ • V^(f ) (101) 

where ^ = ^^70 with the 4x4 matrix 70 = 0-3 (g) as (in this section we put h = 1). We further have 7 = (71,72) = 
(— 2(72, icTi) (g)cr3. The denote the usual Pauli matrices such that {7^, 7^^} = 2g^^jyl4x4^ 9n,u denoting the Minkowski 
tensor where ^0,0 = 1, 9i,i = — 1 with i = 1, 2, and zero otherwise. 

The long range Coulomb interaction in terms of the Dirac spinors reads 

ifee = ^ / cfxcfix'^{xho'^{x)—^^{x')jo'^{x') (102) 

47r J \x — 

where g = e'^ /vp is the dimensionless coupling constant. 

In order to describe disorder effects, the Dirac spinors are coupled to a gauge field A(x), 

Hd^sorder = ^ / {x)T^ (x) A{x) (103) 

where vr characterizes the strength and the 4x4 matrix F the type of the vertex. In general, A{x) is a quenched, 
Gaussian variable with the dimensionless variance A, i.e., 

{A{x)) = , {A{x)A{x')) = ArS\x-x') . (104) 

Inclusion of long ranged correlated disorder is a much more difficult task that has been addressed in |103( I177[ 1178] . 

A perturbative renormalization group (RG) approach to the problem as the one performed in |174H176j allows to 
study the phase diagram of the system as a function of the Coulomb interaction parameter and the strength of the 
various disorder couplings. It is known that the (unscreened) Coulomb interactions in graphene induce an upward 
flow of the Fermi velocity that grows at low energies [23l 1179] . That makes the effective coupling ^ ~ ^ to flow to 
zero. Typically the inclusion of disorder changes the flow of the Fermi velocity due to wave function renormalization, 
i.e., Zy^ Zyp./Z^. In the simplest one loop analysis done in jl74j . from the /3- function /^^^ = Kd\ZypVF^ one 
obtains flow equations for the effective Fermi velocity v^J^ of the form: 



d v'J^ _ 1 



di vp 167r 




(105) 



where = In A/Aq ~ l/e. 

The phase diagrams obtained in |174j for the various types of disorder considered can be summarized as follows: 

i) For a random chemical potential (F = 70), = ^1 remains constant under renormalization group transformation. 
There is an unstable flxed line at Vp = v\/S./ (2e^). In the (^, A)-plane, the strong-coupling and the weak-coupling 
phases are separated by a hyperbola, with the critical electron interaction g* = e'^/vp = 2e^ /{vfA). 
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ii) A random gauge potential involves the vertices F = 271, 1^2- The vertex strength renormalizes as vy = vp. There 
is an attractive Luttinger-like fixed point for each disorder correlation strength A given by = 2e^/A or ^* = A/2. 

iii) For a random mass term F = l4x4, topological disorder F = ^75, and F = ^75, we have '^r = Vp/vs. There is 
thus again an attractive Luttinger-like fixed point for each disorder correlation strength A given by v*p = i/2vieyA 
or g* = {/Aey{2vl). 

More sophisticated calculations do not introduce essential changes in this simple description. 



E. Scattering of charge carriers by gauge fields 

High electron mobility found already in the first work on graphene [13 is one of its most attractive features, in 
view of potential applications. Submicron mean free path is routinely achievable; even order of magnitude higher 
mobility /i :^10^ cm^/Vs was reached recently for freely suspended graphene membranes 1180(1181 . Despite numerous 
experimental and theoretical efforts, the physics of this high mobility is not well understood yet and there are still 
controversial views on the main mechanism limiting the electron mobility. Following a general line of our review, we 
discuss here the role of gauge fields created by corrugations of graphene in its electron transport properties. 

We proceed with a simple and physically transparent semiclassical picture based on the Boltzmann kinetic equation. 
Its applicability to graphene is not a priori obvious, due to potential role of quantum relativist ic effects known as the 
Zitterbewegung (in terms of solid state physics - interband processes mixing electrons and holes) jl82H183j . However, 
it has been formally demonstrated in Ref. |183j . using short-range scatterers as an example, that the interband 
scattering processes are negligible, and, thus, the semiclassical Boltzmann equation is justified, assuming that 

SFT/h:^ l/\\n{kFa)\ , (106) 

where Sp and kp are the Fermi energy (counted from the Dirac point) and wave vector, respectively, and r is the 
electron mean free path related with the resistivity p by the Drude formula |184j 

o 1 

(107) 



eHlN{eF)T{kF)' 



The condition (106) excludes only a relatively small doping interval where the conductivity is close to a minimal 
metallic conductivity of order of /h. The expression for r depends on the scattering mechanism. For point defects 
with the concentration riimp and angular dependent scattering cross section a (0) one has |184[ 1185] 

= riimpVpcrtr (^) , (108) 

where 

C7tr = J dO^^ {1- COS 0), (109) 


is the transport cross section. 

For the case of radially symmetric potential the scattering cross section can be expressed in terms of the scattering 
phases (5^ as [271 [28l [ESI fT86HT90) 

4 

Cr^^ = - ^ Sin^ {Sm - Sm^i) . (110) 

Consider first the case of short-range scatterers with the radius of potential R much smaller than the electron de 
Broglie wavelength = 27: /kp- Then only 5-scattering survives (m = 0), with 5q (k) ~ kR in a generic case, with a 
negligible contribution to the resistivity. 



P - -^UimpR^ (111) 
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(we keep here and further a factor 4 to remind that there are four channels of conductivity, due to two valleys and 
two spin projections). This is not surprising: the massless Dirac fermions in graphene has the same dispersion law as 
light and, thus, obstacles with sizes smaller than the wavelength are inefficient, like in optics |191j . On the contrary, 
for conventional nonrelativistic two-dimensional electrons even weak scattering leads to a formation of shallow bound 
states and to singularities in the scattering matrix [1 92 , namely, Sq {k) ^ 1 / In (kR) . For the case of graphene, such 
situation takes place for a special case of "unitary" , or "resonant" scatterers with a quasibound state located near the 



Dirac point [ 27tll69[ll86j . In such a case, the contribution to the resistivity is much stronger than (111) |27lll93[ll94| . 



where n is the charge carrier concentration. 

For the case of long-range Coulomb scattering potential, the scattering phases are energy independent jl87H189|ll95j 
which leads to a contribution to the resistivity inversely proportional to n; taking into account the screening does not 
change this result |177( I196| 1197] . A rough estimation for the resistivity reads [ 197j 

p^2o4^ (113) 

Experimentally [14l [15], the resistivity is roughly inversely proportional to the charge carrier concentration, or, 
equivalently, the mobility is independent, or weakly dependent on n. Since the scattering by Coulomb charges is 
the simplest mechanism explaining this behavior it is not surprising that this was considered as an explanation 
"by default". This explanation seems to be in agreement with experimental data on graphene chemically doped by 
potassium |198j . On the contrary, experiments on the chemical doping by gaseous impurities such as NO2 show a 
rather weak dependence of the mobility on charge carrier concentration [199^. Recently, it was demonstrated that the 
mobility is also weakly sensitive to dielectric constant of substrate or to coverage of graphene by polar liquids with 
high such as ethanol or water [200] . Clusterization of the charge impurities due to low diffusion barriers may be a 
possible explanation of their relatively small contribution to the resistivity ^201(1202] . Anyway, the experimental data 
[200 probably mean that charge impurities cannot be the main factor limiting electron mobility in graphene. 

Resonant scatterers such as vacancies [TW or covalently bond impurities [202", 203^ might be another option. 
However, vacancies in graphene have huge energy, about 7 eV, and should be rather exotic (if not created intentionally 
as in experiments |204j ): as for the covalent impurities it is not clear why they should always have resonant levels 
close to the Dirac point. For these reasons, an alternative scenario [194 , the scattering by gauge fields created by 
frozen ripples, is worth to be considered seriously. 

Let us consider a generic perturbation of the form 

pp' 

where is the Dirac spinor dependent on quasimomentum p and 

contains both scalar (electrostatic) and vector (gauge field) random potentials. Following Ref. [194j we will use Born 
approximation in the perturbation H' and Kubo-Nakano-Mori approach |205ti207] (as a recent example illustrating 
technical details see, e.g., Ref. [208[). This approach allows us to obtain, in a simple and straightforward way, the 
results which are equivalent to solution of Boltzmann equation by variational principle [184^ . The inverse mean free 
path time is given by 



1 1 



01) 



00 

J dt{[j,{t),H'{t)][H',j,]), (116) 



where jx is the current operator in x direction. Use of this approach in graphene should be done carefully, namely, 
only intraband (electron-electron or hole-hole, depending on the positions of the Fermi energy) contributions to the 
operators jx and H' should be taken into account. An accurate justification of such procedure is done in Ref. |183j . 
For the case of static disorder, the result takes the form 

1 47r — ^ 2 2 

hN(e ) Z^^(^P ~^-^)^(^P^ -gF)(cos(9p-cos(9p/) |iypp/| , (117) 
pp' 
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where 9p is the polar angle of the vector p and 



pp' 



1 



(0) 



1 + exp [- 



^p')] 



{4t + iA^^i,) exp (-i^p) + - i^SJ,) exp (ie^,) 



(118) 



Note that the electrostatic part of the scattering disappears for the back scattering (^p — ^p/ = tt) which is related to 
the "Klein tunneling" [25 but this is not the case for the scattering by vector potential. For a rough estimations in 
order of magnitude Eq. (117) can be rewritten as 



27r7V {ep) 

h 



v_ 



{AqA_ 



(119) 



where q = p'— p is the scattering vector. Sin ce N {sf) cx kp^ to have concentration independent mobility one needs 
that one of the correlation functions in Eq. (119) scales as Note that the second term, that is, the scattering 



by gauge fields, is not sensitive to dielectric screening. 

Cons ider n ow the case of scattering by gauge fields created by intrinsic ripples due to thermal fiuctuations (see 
Section IV D). For not too small wave vectors, in the harmonic regime, q > one can neglect a coupling between 



bending and stretching modes. This case is relevant for electron transport assuming that 



(120) 



which does not look too restrictive. Further we will assume this condition to be satisfied. Thus, the vector potential 
is quadratic in derivatives dh/dx^dh/dy and the estimation of the correlation function in (119) is 



(AqA_ 



/ hVF 

V a 



q-qi^qi^-q+q^/i-qs) [(q " qi ) ' qi] [(q " q2) • q2] , 



(121) 



and the correlator in right-hand side of Eq. (121) can be decoupled by Wick theorem since the field h (x^y) can be 



considered as a Gaussian in this regime. Thus, the scattering rate under the condition (120) is determined by the 
correlation function (/iq/i-q) , that is, a Fourier component of the correlator 



r{r) = {[h{r)-h{0)f) 



For thermally excited ripples in the harmonic regime 

r(r) 



(122) 



(123) 



which is a Fourier transformation of the Eq. ( |82| ). However, we consider first a more general case F (r) ex r^^ (for 
instance, for the ripples due to roughness of substrate one could expect H ^ 1/2, see Ref. [45^). 
For 2H < 1, the correlation function (121) has a finite limit at g = 



(AqA_q)^^o« i a) i?2' 



(124) 



where z and R are the characteristic height and radius of the ripples, respectively. The corresponding contribution 
to the resistivity is 



h 



4e2 i?2a2 ■ 



(125) 



For 2 > 2H > 1, the resistivity is proportional to n^~^^ and, for 2H = 1, to |ln(/ci?a)| . In all these cases the 

corresponding contributions are too small in comparison with experi menta lly observed resistivity of graphene^ 

Let us now come back to the case of thermally excited ripples (123), H = 1. The integral in Eq. ( |121[ ) is 
logarithmically divergent at small wave vectors, but this divergence should be cut at qi q* . The corresponding 
contribution to the resistivity is estimated as 



h 

4^ 



TV \\n{q*a) 
na J n 



(126) 
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This contribution gives us correct concentration dependence of the resistivity and, for room temperature T = 300K, 
correct order of magnitude of mobihty, /i o:^ 10^ cm^/Vs. However, it predicts a strong temperature dependence of 
the resistivity, whereas experimentahy it is very weak for the case of graphene at substrate |2Q9j . The authors of Ref. 
[194 postulated that the ripple structure for graphene on substrate is frozen and becomes temperature independent 
below the room temperature. 

The status of this hypothesis is not clear yet, neither theoretically nor experimentally. Whereas earlier STM studies 
of the ripples of graphene on substrate [^S*, "46 have found corrugations more or less similar to those of the substrate, 
recent work |21Qj claims that an "intrinsic" component postulated in Ref. [194 is also noticeable. Further investigations 
of the temperature dependence of the ripple structure are desirable to clarify the issue. From the theoretical point of 
view, it was demonstrated in Ref. [155 that ripples can attract adatoms and chemical groups (such as hydrogen or 
hydroxyl) and can be stabilized by them. Again, the problem requires more studies to clarify the situation. 

Up to now, we consider out-of-plane deformations as a classical static field. In terms of quantum mechanics, the 
corresponding scattering mechanism is described as the two-phonon scattering processes by fiexural phonons ^209j. 
Introducing in a standard way |184j phonon creation and annihilation operators 6^, 6q one can represent the out-of- 
plane displacement field as 



h 



2MuJr 



(127) 



where M is the mass of carbon atoms, Wq = ^/n/pmO^ is the frequency of flexural phonons, Pm is the mass density. 
Since the vector potential is quadratic in dh/dx,dh/dy the interaction Hamiltonian has the form 



^' = Y1 -^Pqkcj,Cp+k_q (5Lq + 6q) (bl + 6_k) , 



(128) 



pqk 



where we take into account only intraband scattering processes (cj,, Cp are electron creation and annihilation operators) 
and the amplitude A behaves at small wave vectors as 



^pqk 



qk. 



(129) 



Then, we simply substitute the Hamiltonian (128) into Eq. (109). The calculations are very similar to those for two- 



magnon scattering processes in half-metallic ferromagnets [2Q8j . The result (where we skip for simplicity numerical 
factors of order one) reads 
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(130) 



where fp = f (Sp) and A^q = Nb (c^q) are the Fermi and Bose distribution functions, respectively. For temperatures 
T > hw2k —IK phonons can be considered as classical neglecting the phonon frequencies in the conservation laws 



and replacing N^^ by T/huOd in Eq. (130). Then, the result coincides with Eqs. (119), (126). The case of very low 



temperatures was considered in Ref. [211^. The minimal conductivity from a geometric model of topological lattice 
defects has been computed in |103j . 



VI. OBSERVABLE CONSEQUENCES OF THE GAUGE FIELDS: MESOSCOPIC EFFECTS, STRAINS 

IN SUSPENDED SAMPLES 



A. Aharanov-Bohm phases in suspended samples 



Mesoscopic deformations have been observed in graphene, and they can be induced in a controlled way. A par- 
ticularly appropriate setup is a suspended graphene flake [43l |45l 11801 11811 l212H219j . These systems are being very 
actively investigated, as the carrier mobility can be substantially larger than that of graphene on a substrate. The 
conductance minimum at the charge neutrality point is also sharper than in non-suspended samples. It has been 
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observed that these samples are under tension, which can be modulated by the electric force between the flake and 
the gate. 

Typical sizes, L, of suspended samples are of the order of a few microns. Away from very close to the neutrality 
point these dimensions are much larger than the Fermi wavelength and the screening length, which is proportional to 
it. Hence, the long wavelength scalar potentials, Vo{f), which may be generated by compression [62, 63 are screened. 
A simple Thomas-Fermi approximation leads to Vscr{r) ~ Vo{r)/{kFTL) where kpT = e^kp / {tivf) is the screening 
length. Gauge flelds, on the other hand, remain unscreened. 



1. Electronic transport in ballistic samples 

We flrst assume that the samples are free of disorder and the only source of electron scattering is the gauge fields 
induced by the deformations. The simplest geometry which can be studied is a rectangular flake, as sketched in Fig. 
[t] If the bending energy is neglected, a complete analytical form can be found for the deformation [121]: 

T / N ho f 2 



L2 V 4 

where E = 4(A + /i)/i/(A + 2/i) ^ 24 eV A~^ is the Young's modulus of graphene. The in-plane stresses are: 

+ i (132) 

where the boundary condition Ux{^L/2) = has been used. As a result, the only non zero component of the strain 
tensor is: 

= ^ (133) 

The effect of the uniform strain is to induce a constant gauge field in the suspended region. Ay = Puxx/d- The 
effect of this field is to shift the momentum parallel to the interface between the suspended and clamped region, ky. 



An electron outside the suspended region with momentum k = {kx^ ky) and energy eg = vpyk'^ + ky can be either 

refiected at the boundary, with momentum k' = {—kx^ ky) or transmitted, with momentum k" = (/c^, ky -\-Ay) with k'^ 
such that /c^ {ky -\- Ay)'^ = k^ -\- ky . For incident momenta in the range k^ < 2Ayky -h Ay there are no transmitted 
states, and the incident electron is refiected. The transmission amplitude, T, is determined by the equations: 

1-R = T 

kx -\- iky —kx-\-iky ki^ -\- i{ky -\- Ay) /iq/ia 

+ — k — ^ = k ^ 

and the transmission coefficient is T = \T\'^{k'^/kx). At normal incidence, ky = 0, the transmission coefficient can be 
expanded, T ^ 1 - Al/{4kl). The fuh expression for the transmission coefficient in the setup shown in Fig. 7 with 
two boundaries, can be obtained analytically [121^. The presence of two barriers leads to Fabry- Perot interferences 
as function of the incident angle, = 8iTctdin{ky / kx) . The angular dependence of the transmission is shown in Fig. [s] 
Note the skew scattering in the two valleys. 

The above analysis is only valid for clean, ballistic systems. A finite concentration of impurities breaks the conserva- 



tion of parallel momentum, which leads to the suppression of the transmission through each interface, Eq. (134). The 



lack of momentum conservation gives a finite transmission when strains in a ballistic sample completely suppresses 
it [121\ Evanescent waves allow carriers to penetrate in the suspended region. These carriers can be scattered into 
propagating modes by the impurities. A perturbative calculation shows that the induced transmission grows linearly 
with the concentration of impurities. The effect is proportional to the scattering strength, which can be parametrized 
in terms of the mean free path induced by the impurities. 



B. Interferences between strains and real magnetic fields 



The gauge fields associated to strains interfere with real magnetic fields applied to the system. When the strain can 
be treated as a small perturbation, the resulting effect is the existence of an inhomogeneous effective field distribution. 
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FIG. 17: Sketch of the expected effect in the presence of strains and a constant magnetic field. The geometry is the same as 
in Fig. [7]. The strains induce currents which lead to backscattering between the chiral edge states. A small concentration of 
impurities suppresses the effect. 



given by the sum of the constant applied field and the effective field induced by the strains. Electronic states in the 
vicinity of the inhomogeneities are deformed, and can lead to localized "snake states" , where current fiows along the 
perturbed region, as in the presence of a nonuniform magnetic field |22QH222 j. 

The qualitative effect of strains in the geometry in Fig. [7] is shown in Fig. 17 The strains induce currents parallel 
to the interface between the suspended and the non suspended regions, which allow for backscattering between the 
chiral edge states created by the magnetic field [223]. The states along the interface are not chiral, and impurities 
induce backscattering between them. If this effect is sufficiently strong, these states will become localized, and the 
coupling between the chiral edge channels will be suppressed. 

The hamiltonian which describes the system, in the absence of impurities, and neglecting the edges, is: 



1-L = vf 



id:, - iky ± iA^^^^^^6>(x) 



- iky T iAl^'-'^'^'e^x) 



(135) 



, where Ib is the magnetic length associated to the external fiel d, an d we assume that, in a system infinite in the y 
direction the momentum ky is conserved. The two signs in Eq. (135) refer to the two valleys in the Brillouin Zone. 
By squaring the hamiltonian, one obtains the effective Schrodinger equation: 



.,2 q2 



strain 
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(136) 



where X- = kyl'^ and x+ = x_ + A^^^^*^. For a given value of ky, the eigenvalue Cky depends only on vp and 
Ay^^^'^^ls' The strain induces two effects: i) it shifts the centers of the Landau levels, and ii) a delta function potential 
appears at the interface. At distances from the interface greater than Ib, the Landau levels are not perturbed, and 
:^VF^/n/lB-, where n is the Landau level index. For Ay^^'^'^'^lB <C 1 we can use perturbation theory to analyze 
the shift in the energies of the Landau levels. The leading effect is due to the delta function potential. The maximum 
change in the energy is: 



it 



VFyn 
Ib 



l + ±c^ 



Ar,,l 



2n 



(137) 



where is a numerical constant. The Landau levels become dispersive, with a velocity 'U* = de^^/dky ~ 

-^vpA^y^^^^lB- The Landau level energy has a maximum, or a minimum, for a value of ky such that the center 
is close to the interface. Hence, there are two counter propagating states per Landau level and per valley. These 
modes have a width ^ Ib- Using the Born approximation, it can be shown that the reffection coefficient induced by 



a single charged impurity is of order e /vp for A 



strain 



1. 
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FIG. 18: Sketch of applied strains which lead to an almost constant effective magnetic field inside the sample. 



C. Generation of effective magnetic fields. Strain engineering 



The previous analysis suggests that strains can generate a wide variety of effective magnetic fields. This has 
open the possibility of using "strain engineering" in order to modify the electronic properties of graphene samples 
[84 1 \112\ I155[ l224H233j . Changes in the strain can enhance the chemical activity of graphene, or confine the massless 
carriers, avoiding the Klein paradox [25]. Experimentally, strains can be induced in graphene flakes by bending 
the substrate |234H236j . and the amount of strain can be measured by Raman spectroscopy. Strains can be due 
to the difference in expansion coefficients between graphene and the substrate, leading to wrinkles [217] . There is 
experimental evidence that the electronic properties are modified in strained samples [237] . 

A particularly interesting possibility is the generation of an effective constant magnetic field. The results in the 
Appendix [B] show that a strain distribution with trigonal symmetry induces such a field. An experimental setup which 
leads to an approximate uniform field is shown in Fig. 18, see [227[ 1232] . Such situation will change drastically the 
electronic properties of the system. The electronic spectrum will be split into Landau levels, with counterpropagating 
chiral modes at the edges. The combination of insulating behavior in he bulk and gapless modes at the edges is 
reminiscent of a topological insulator [22l 1238] . 



D. Effective electric fields 



A time-independent strain gives rise to static gauge fields. Hence, a time-dependent strain will induce effective 
electric fields [239 . These fields appear in oscillating graphene sheets [ISHl l24Qj . As in the static case, the fields 
acting on different valleys have opposite sign, so that there is no net charge transport and the field is not screened. 
The time dependence of the scalar potential also induced by the strains [62] also leads to an effective electric field, 
although with the same sign for the two valleys. This potential can be screened. The effective electric fields induce 
electronic currents, which dissipate energy. This ohmic damping reduces the quality factor of the oscillator. 

The description of the decay of long wavelength phonons into electron-hole pairs by means of an effective electric 
field acting on the electrons has also been used to understand phonon damping in carbon nanotubes ^241j (for a 
review of the extensive literature on nano electromechanical systems, NEMS, based on carbon nanotubes, see |242j ). 
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In a clean system, the losses due to currents induced by the effective electric fields are equivalent to the phonon 
lifetime which can be calculated using second order perturbation theory. A formulation using effective fields can 
easily be extended to diffusive systems, where the electronic wave functions can only be estimated using semiclassical 
arguments. Standard perturbation theory is adequate for clean, ballistic carbon nanotubes, but a description using a 
finite mean free path and local conductivities is more appropriate for large graphene samples, as the elastic mean free 
path is typically much smaller than the dimensions of the sample. It is worth noting that the existence of two counter- 
propagating valley currents imply that the ohmic dissipation in a graphene resonator is limited, at high temperatures, 
by valley drag, which, in turn, is dependent on the strength of the electron-electron interaction. 

VII. CONCLUSIONS 

To summarize, we have presented an overview of the role of gauge fields in graphene and other carbon-based 
materials. The geometric structure of honeycomb lattice naturally leads, in the continuum limit, to massless Dirac 
fermions as elementary excitations. The internal degree of freedom which plays the same role as spin in quantum 
electrodynamics has here a very simple geometric meaning, it is just the sublattice label (pseudospin). Additional 
internal quantum numbers are the valley label (isospin) and real spin. If we neglect spin-orbit coupling (which is 
very small in graphene) the latter is split from other degrees of freedom. We distinguish two types of gauge fields: 
topological, associated to lattice defects, such as disclinations and dislocations, and related to modulations of the 
hoppings, induced by elastic deformations. 

The coupling of the topological fields to the electrons does not depend on the parameters of the material. Topological 
fields come in two mutually non commuting varieties, depending on whether they mix valleys or not. The gauge fields 
associated to long wavelength strains do not mix valleys, and their coupling to the electrons depend on the value of 
the strength of the electron-phonon interaction. 

Being exactly a two-dimensional system, graphene shows a thermally induced bending instability which leads to 
intrinsic ripples (corrugations) even in a free sample. For the case of exfoliated graphene on a substrate, the ripples 
will arise also due to interaction with the substrate. In both cases, the ripples are sources of a random pseudomagnetic 
gauge field which modifies the electronic transport in the weak (ant i) localization regime, and it also induces essential 
contributions to the resistivity. The ripples may be one of the main limiting factors restricting mobility of charge 
carriers in graphene. 

Regular gauge fields, or even artificially designed fields can exist in suspended samples, where they can lead to 
a number of novel features unique to graphene, such as the confinement of electrons, or the formation of localized 
orbit als similar to the Landau levels which exist in real magnetic fields. 

The study of gauge fields in graphene establishes important relations with contemporary mathematical physics. 
For example, due to the Atiyah- Singer index theorem the gauge fields in both single-layer and bilayer graphene can 
create states with zero energy, the states being chiral (which means that for a given valley they belong to only one 
sublattice) . The inhomogeneity of pseudomagnetic field created by the ripples should broaden all Landau levels in the 
quantum Hall regime except the zero-energy one. This prediction of the index theorem can be checked experimentally 
and seems to be in an agreement with the available data. It is interesting to note that the gauge potentials which arise 
from lattice deformations have intrinsic meaning, although their effects on the electrons are invariant under gauge 
transformations, so that different physical deformations can lead to the same electronic structure. 

We have not explored here extensions to other systems, such as multilayered graphene, or the topological insulators. 
The same arguments used to define gauge fields in single layer graphene also apply to multilayered samples. The 
effect of a given deformation, or gauge field, on the electrons will depend on the number of layers and in the way they 
are connected. In topological insulators with surface states described by the Dirac equations, gauge fields can also be 
defined. Because of the different symmetries of these compounds, such fields should arise from perturbations which 
break time reversal symmetry. 
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Appendix A: The Dirac equation in curved space 



In general when trying to formulate a mathematical relation defined in a fiat space to a curved space, one uses 
a general covariance principle that amounts to substitute any given magnitude transforming as a tensor in the fiat 
space by the corresponding magnitude transforming as a tensor under general transformations in the curved manifold. 
This simple substitution is complicated for spinors because there are no spinor representations in the group of general 
transformations. This makes necessary to introduce an alternative formalism based on tetrads |95lESl- Instead of 
the usual metric ^^^^ we must introduce at each point X described in arbitrary coordinates, a set of locally inertial 
coordinates and the fielbein fields e^^{x)^ a set of orthonormal vectors labelled by a that fixes the transformation 
between the local (Latin indices) and the general (Greek indices) coordinates: 

We can now list the various geometric objects needed to derive the Dirac physics in the curved space: The metric 
tensor of the curved manifold g^u{x) is related to the fiat, constant metric rjab by the equation 

g^,u{x) = e''^{x)el{x)r]ah] (A2) 
its determinant, needed to define a scalar density lagrangian is given by 

v^=[det(5^,)]i/2 = detK(x)]. (A3) 



The curved space gamma matrices 7^(x) satisfying the commutation relations 

{r{x)Y{x)} = 2g^%x), (A4) 

are related with the constant, fiat space matrices 7^ by 

7^(x) = e^^(x)7". (A5) 

The most complicated object needed to complete the analysis is the spin connection Vt^{x) associated to the spinor 
covariant derivative. This is an important object for the physics since it acts as a gauge field. The construction of 
the spin connection is done by observing that the derivative of the spinor does not transform as a vector under a 
coordinate transformation in the fiat tangent space. It does if we introduce for the spinors a covariant derivative of 
the form 



d 

dx^ ^ ^ 



(A6) 



which has a standard transformation under a change of coordinates in both the vector and the spinor indices. To get 
the structure of the l^^(x) matrices it suffices to consider their transformation properties in ffat space. For the spin 
one half representation they take the form 



1 

4 

with 



= -jalbel{x)g^^x)VM^), (A7) 



where P-^^ is the usual affine connection which is related to the metric tensor by 



V^e^ = d^el - r'el (A8) 



'"^ " 2^ ^dx^^ ^ dx- dx- ^ ' 

Finally the Dirac equation is given by 

i-i^{x)[d^ + ^^{x)]^{x) = 0. (AlO) 

We will exemplify the formalism by giving the details of the computation of the various geometrical factors for the 
two main problems described in Section [llH the topological defects and the smooth gaussian deformation. 



38 



We begin by describing the formalism for a smooth deformation where the curvature is non singular everywhere. 
We start by embedding a two-dimensional surface with polar symmetry (this is only for simplicity and can be easily 
extended to any shape) in three-dimensional space (described in cylindrical coordinates). The surface is defined by a 
function z{r) giving the height with respect to the fiat surface z=0, and parametrized by the polar coordinates of its 
projection onto the z=0 plane. The metric for this surface is obtained as follows: We compute 



dz^ = {^fdr^ = af{r)dr^ 

and substitute for the line element: 

ds^ = dr^ + r'^de^ + dz"^ = (1 + af{r)) dr^ + r'^de^ 
In particular the the gaussian bump shown in Fig. |6]is defined by: 

z = Aexp(-r^/6^), 

so that 

dz^ = ^4:r^exp{-2r^/b^)dr^, 



(All) 

(A12) 
(A13) 

(A14) 



which corresponds to Eq. (A12) with 

a = (A/b)^ , /(r) = 4(r/6)2exp(-2rV^'). 
From the line element we can write the metric in a more usual form: 

-(l + a/(r)) 







(A15) 



where we have omitted the time coordinate which is tri vial in the case of a static background. The affine connection 
which only depends on the metric is for the metric (A15): 
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where = df /dr^ and the rest of the elements are zero or related by symmetry. 
The geometrical (gaussian) curvature K is 



K{r) = 



af{r) 



2r(l + a/(r))2 



The fielbein fields e^^ satisfy: 



9nv — ^ /i^ vVah^ 



(A17) 



(A18) 



where 77^5 is the identity matrix in two dimensions (note that this relation does not fix e^^ uniquely). We choose the 



e^,, to be 



e\ = {1 -\- afy^'^ cosO e\ = —rsin 



e^Q = r cos 6; 



that reduce to the fiat set when a = 0. Now we can compute the spin connection coefiicients, 



which are found to be: 



12 



l-{l + af)-'/\ 



(A19) 
(A20) 
(A21) 
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the rest being zero or related by symmetry (the spin connection uj is antisymmetric in the upper indices j244j ). 
The spin connection 



% = \^tbaM. (A22) 
turns out to be 

= 0, ne= (A23) 



In the case of topological defects, for the metric defined by (44): 



9^lv 



-1 











g-2A 











g-2A 



(A24) 



the gamma matrices and the spinor connection in the curved background are immediately found to be 



7°(r) = 7° , 7'(r) 



=A(r) 



y (i = i,2) 



and the determinant of the metric tensor is 



g = e 



-2A 



The special feature of the conical defects is that the Riemann curvature tensor is zero everywhere but for the apex of 
the cone where it has a delta function singularity. 

Appendix B: In plane strains and effective fields in radial coordinates 

We analyze a circular graphene quantum dot. We first analyze the distribution of stresses. We assume that the dot 
is free of lattice defects. The stresses are due to forces applied at the boundaries. 
The elastic free energy in circular coordinates is: 



T = — I 27rrdrd0 ( drUr 

2 7 V r 



I 

ji j 27TrdrdO 
drUe 



r J 



idrUr) + \ 

V r r 



(Bl) 



where Ur and uq are the two components of the displacement vector in circular coordinates. The stress tensor is: 



Grr = A ( drUr + ^ H ^ ) + 2/i (drUr) 



aee = drUr + 



deue , Ur 



2/i 



dpue ^ Ur_ 
r r 



CFrO = 2/i 



^eur , ^ ue 

h OrUo 

r r 



(B2) 
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ae/ttr 



B{r,e) 



-i{m — 2)0^m — 2 



• A(m+l)+/x(m+3) 
^ A(m — l)+/j,(m — 3) 
• A(m+l) + /x(m+3y 



Aim(m — l)e 



-m + 1 



4im(m-l)(m-2) ,^^_,^+;:^^_3^ e 



' A(m-l) + /x(m-3) 



-4im{m — l)(m — 2) 



A+/^ 

A(m-l)+/x(m-3) ^ 



i{m — 2)e^m — 2 



m+ 1 



-4im{m — l)e 



TABLE I: Displacements and induc ed e ffective magnetic fields for solutions of the equations of elasticity with well defined 
symmetries under rotations, see Eq. (B4). 



The resulting equilibrium equations for Ur and uq are: 



+ 2/i 
= 



A 



drdeUr doUr dgUe 



2/i 



= 



3deUr dnUe drUo uq drdgUr d^ue 



2r2 



2r 2r2 2r 



We look for solutions of the type: 



Inserting these expressions in Eqs. |B3[ we find: 



-m 



Ug 



(B3) 



(B4) 



l) + /i (-2m^ + 2 + n^)] a^- 



+(m) [A(-m + 1) + /i(3 - m)] = 
(— m) [A(m + 1) + /i(3 + m)] 

[An^ + /i(2n^ - + 1)] = 



(B5) 



These equations admit solutions for n = ±(m ± 1). For m = there are only two solutions, with n = ±1. For m = 0, 
a finite value of ag leads to a rotation of the disk, and does not change the physical properties. 

The induced gauge field on the Dirac electrons of the graphene layer depends on the relative orientation of the 
graphene axes. Choosing the direction 6> = as one of the lattice axis, the gauge potential can be written as: 



Urf 

r r 
dpue _ Ur 
r r 



cos(3i9) - 
sin(3(9) 



9Ur 



r 

dpUr 

r 



drUo - 

+ drUe 



ue_ 
r 



ue_ 
r 



sin(3l9) 
cos(3i9) 



(B6) 



where /? = d\og{t) / d\og{a) « 2 — 3 is the logarithmic derivative of the nearest neighbor hopping t with respect to the 
nearest neighbor distance, o. The effective magnetic field is: 



B 



deAr 



drAg 



Ae 



(B7) 



The effective magnetic fields induced by solutions of equilibrium elasticity of the type described in Eq. ( B4 ) are given 
in Table [l]. Non trivial solutions which induce a constant field can be obtained by setting m = 2 and combining the 
first and second solutions in Table [l]. 
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